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I. INTRODUCTION 

Dynamics of standard and compact astrophysical bod- 
ies could be aimed to unravelling the information con- 
tained in the gravitational wave (GW) signals. Several 
situations are possible such as GWs emitted by coalescing 
compact binaries-systems of neutron stars (NS), black 
holes (BH) driven into coalescence by emission of grav- 
itational radiation (considered in the inspiralling, merg- 
ing and ring-down phases, respectively), hard stellar en- 
counters, and other high-energy phenomena where GW 
emission is expected. Furthermore the signature of GWs 
can be always determined by the relative motion of the 
sources. In this review paper, we want to discuss the 
problem of how the waveform and the emission of GWs 
depend on the relative motions in Newtonian, Relativistic 
and post-Relativistic regimes as for example situations 
where gravitomagnetic corrections have to be considered 
in the orbital motion. 

As a first consideration, we have to say that the prob- 
lem of motion, i. e. the problem of describing the dynam- 
ics of gravitationally interacting bodies, is the cardinal 
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problem of any theory of gravity. From the publication 
of Newton's Principia to the beginning of the twenti- 
eth century, this problem has been thoroughly investi- 
gated within the framework of Newton's dynamics. This 
approach led to the formulation of many concepts and 
theoretical tools which have been applied to other fields 
of physics. As a consequence, the relationship between 
Einstein's and Newton's theories of gravity has been, and 
still is, very peculiar. On the one hand, from a theoretical 
point of view, the existence of Newton's theory facilitated 
the early development of Einstein's theory by suggesting 
an approximation method (called post- Newtonian (PN)) 
which allowed to draw very soon some observational con- 
sequences of General Relativity (GR). Indeed, the PN ap- 
proximation method, developed by Einstein himself [ij, 
Droste and de Sitter 0, Q within one year after the pub- 
lication of GR, led to the predictions of i) the relativistic 
advance of perihelion of planets, ii) the gravitational red- 
shift, in) the deflection of light, iv) the relativistic preces- 
sion of the Moon orbit, that are the so-called "classical" 
tests of GR. 

On the other hand, as emphasized by Eisenstaedt [1], 
the use of PN approximation method has had, from a 
conceptual point of view, the adverse side-effect of in- 
troducing implicitly a 'neo-Newtonian' interpretation of 
GR. Indeed, technically this approximation method fol- 
lows the Newtonian way of tackling gravitational prob- 
lems as closely as possible. But this technical reduction 
of Einstein's theory into the Procrustean bed of New- 
ton's theory surreptitiously entails a corresponding con- 
ceptual reduction; the Einstenian problem of motion is 
conceived within the Newtonian framework of an "ab- 
solute" coordinate space and an "absolute" coordinate 
time. However, some recent developments oblige us to 
reconsider the problem of motion within Einstein's the- 
ory. On the other hand, the discovery of the binary pulsar 
PSR 1913-1-16 by Hulse and Taylor in 1974 j^, and its 
continuous observation by Taylor and coworkers (see ref- 
erences [1,01) J led to an impressively accurate tracking 
of the orbital motion of a NS in a binary system. This 
means that it is worth reconsidering in detail, i.e. at its 
foundation, the problem of motion also in relation to the 
problem of generation and detection of GWs. In other 
words, the motion of sources could give further signatures 
to GWs and then it has to be carefully reconsidered. 

The first part of this review paper is devoted to the 
theory of orbits. The most natural way to undertake 
this task is starting with the discussion of the Newtonian 
problem of motion then we consider the relativistic prob- 
lem of motion, in particular the PN approximation and 
the further gravitomagnetic corrections. 

The theory of orbits can be connected to GWs since 
studies of binary systems prove, beyond reasonable 
doubts, that such a form of radiation has to exist. Detect- 
ing the waves directly and exploiting them could result 
a very impressive way to study astrophysical objects. In 
other words, the detection of GWs could give rise to the 
so-called Gravitational Astronomy. 



In view of this achievement, it is relevant to stress 
that GW science has entered a new era. Experimen- 
tallyB, several ground-based laser-interferometer GW de- 
tectors (10-1 kHz) have been built in the United States 
(LIGO) 0, Europe (VIRGO and GEO) [H, [HI and 
Japan (TAMA) 12], and are now taking data at designed 
sensitivity. 

A laser interferometer space antenna (LISA) [l^l 
(10~'*-10~^ Hz) might fly within the next decade. 

From a theoretical point of view, last years have wit- 
nessed numerous major advances. Concerning the most 
promising GW sources for ground-based and space-based 
detectors, i.e. binary systems composed of NS, BHs, our 
understanding of the relativistic two-body problem, and 
the consequent GW-generation problem, has improved 
significantly. 

Knowledge has also progressed on the problem of mo- 
tion of a point particle in curved spacetime when the 
emission of GWs is taken into account (non-geodesic mo- 
tion) [3, [ill- Solving this problem is of considerable 
importance for predicting very accurate waveforms emit- 
ted by extreme mass-ratio binaries, which are among the 
most promising sources for LISA |16| . 

The GW community working at the interface be- 
tween the theory and the experiment has provided tem- 
plates fir'Jl^ for binaries and developed robust algo- 
rithms (2Q|, [21)] for pulsars and other GW-sources observ- 
able with ground-based and space-based interferometers. 
The joined work of data analysts and experimentalists 
has established astrophysically significant upper limits 
for several GW sources [22h24| and is now eagerly wait- 
ing for the first detection. 

In this research framework, searching for criteria to 
classify how sources collide and interact is of fundamental 
importance. A first rough criterion can be the classifica- 
tion of stellar encounters in collisional, as in the glob ular 
clusters, and in collisionless as in the galaxies |25| . A 
fundamental parameter is the richness and the density of 
the stellar system and then, obviously, we expect a large 
production of GWs in rich and dense systems. 

Systems with these features are the globular clusters 
and the galaxy centers. In particular, one can take into 
account the stars (early-type and late-type) which are 
around our Galactic Center, e.g. Sagittarius A* {SgrA*) 
which could be very interesting targets for the above 
mentioned ground-based and space-based detectors. 

In recent years, detailed information has been achieved 
for kinematics and dynamics of stars moving in the 
gravitational field of such a central object. The sta- 
tistical properties of spatial and kinematical distribu- 
tions are of particular interest (see e.g. [27l - l29| ). Us- 
ing them, it is possible to give a quite accurate estimate 
of the mass and the size of the central object: we have 
(2.61 ± 0.76) X lO^Af0 concentrated within a radius of 



^ GW experiments started with the pioneering work of Joseph We- 
ber at Maryland in the 60s 
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0.016pc (about 30 light-days) [30|, |31| . More precisely, in 
[30| , it is described a campaign of observations where ve- 
locity measurements in the central arcsec^ are extremely 
accurate. Then from this bulk of data, considering a field 
of resolved stars whose proper motions are accurately 
known, one can classify orbital motions and deduce, in 
principle, the rate of production of GWs according to the 
different types of orbits. 

These issues motivate this review paper in which, by a 
classification of orbits in accordance with the conditions 
of motion, we want to calculate the GW luminosity for 
different types of stellar encounters and orbits (see also 

MM)- 

Following the method outlined in [s^, [1^, we in- 
vestigate the GW emission by binary systems in the 
quadrupole approximation considering bounded (circu- 
lar or elliptical) and unbounded (parabolic or hyperbolic) 
orbits. Obviously, the main parameter is the approach- 
ing energy of the stars in the system (see also [36| and 
references therein). We expect that gravitational waves 
are emitted with a "peculiar" signature related to the 
encounter- type: such a signature has to be a "burst" 
wave-form with a maximum in correspondence of the 
periastron distance. The problem of bremsstrahlung-like 
gravitational wave emission has been studied in detail by 
Kovacs and Thorne [37[ by considering stars interacting 
on unbounded and bounded orbits. In this review paper, 
we face this problem discussing in detail the dynamics 
of such a phenomenon which could greatly improve the 
statistics of possible GW sources. For further details see 

also [H, [3113,13® Hi- 

The review is organized as follows. In Part I, as we 
said, we discuss the theory of orbits. In Sec. 2, we start 
with the Newtonian theory of orbits and discuss the main 
features of stellar encounters by classifying trajectories. 
Sec. 3 is devoted to orbits with relativistic corrections. 
A method for solving the equations of motion of binary 
systems at the first PN-approximation is reviewed. The 
aim is to express the solution in a quasi-Newtonian form. 
In the Sec. 4, we study higher order relativistic correc- 
tions to the orbital motion considering gravitomagnetic 
effects. We discuss in details how such corrections come 
out by taking into account "magnetic" components in 
the weak field limit of gravitational field. Finally, the or- 
bital structure and the stability conditions are discussed 
giving numerical examples. Beside the standard perias- 
tron corrections of GR, a new nutation effect have to be 
considered thanks to corrections. The transition to 
a chaotic behavior strictly depends on the initial condi- 
tions. The orbital phase space portrait is discussed. 

Part II is devoted to the production and signature of 
gravitational waves. We start, in Sec. 5, by deriving the 
wave equation in linearized gravity and discuss the gauge 
properties of GWs. Sec. 6 is devoted to the problems of 
generation, emission and interaction of GWs with a de- 
tector. In Sect. 7, we discuss the problem of GW lumi- 
nosity and emission from binary systems moving along 
Newtonian orbits. The quadrupole approximation is as- 



sumed. Sect. 8 is devoted to the same problem taking 
into account relativistic motion. In Sect. 9, also grav- 
itomagnetic effects on the orbits and the emission are 
considered. In Sect. 10, as an outstanding application of 
the above considerations, we derive the expected rate of 
events from the Galactic Center. Due to the peculiar 
structure of this stellar system, it can be considered a 
privileged target from where GWs could be detected and 
classified. Discussion, concluding remarks and perspec- 
tives are given in Sect. 11. 



Part I 

Theory of orbits 

II. NEWTONIAN ORBITS 

We want to describe, as accurately as possible, the 
dynamics of a system of two bodies, gravitationally in- 
teracting, each one having finite dimensions. Each body 
exerts a conservative, central force on the other and no 
other external forces are considered assuming the sys- 
tem as isolated from the rest of the universe. Then, we 
first take into account the non-relativistic theory of or- 
bits since stellar systems, also if at high densities and 
constituted by compact objects, can be usually assumed 
in Newtonian regime. In most cases, the real situation is 
more complicated. Nevertheless, in all cases, it is an ex- 
cellent starting approximation to treat the two bodies of 
interest as being isolated from outside interactions. We 
give here a self-contained summary of the well-known or- 
bital types [1^, which will be extremely useful for the 
further discussion. 



A. Equations of motion and conservation laws 

Newton's equations of motion for two particles of 
masses mi and m2, located at ri and r2, respectively, 
and interacting by gravitational attraction are, in the 
absence of external forces. 



dpi _ ^ mim2 
dt ^ 
dp2 



dt 



-G- 



mi 7712 



3(i"i -r2), 

-(ri-r2), (2.1) 



dri 

where pi = momentum of particle i, {i 

1, 2), and G is Newtonian gravitational constant. 

d 



dt 



(P1+P2) = 0, 



or, with P = Pi + P2 denoting the total momentum of 
the two body system, 

P = const . 
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Thus we have found a first conservation law, namely 
the conservation of the total momentum of a two-body 
system in the absence of external forces. We can make 
use of this by carrying out a Galilei transformation to 
another inertial frame in which the total momentum is 
equal to zero. Indeed, let us apply the transformation 



hence pi 



Pi 



and if we choose v 



= ri-vt, i = l,2 
m-iV and hence, with M — nii 

> P' = P - Afv , 
P 



M 



, then the total momentum is 



equal to zero in the primed frame. We also note that the 
gravitational force is invariant under the Galilei transfor- 
mation, since it depends only on the difference ri — r2. 
Thus let us from now on work in the primed frame, but 
drop the primes for convenience of notation. We can now 
replace the original equations of motion with the equiv- 
alent ones, 

P = 0, ^ = -G^r, 
at r-^ 

where r = ri — r2, r = |r|, and p = Pi — P2- Next we 
introduce the position vector R of the center of mass of 
the system: 



R 



miTi + m2V2 

mi + 7712 



hence 



P = M- 



~dt 



and hence from P we have 

R = const . 

and we can carry out a translation of the origin of our 
coordinate frame such that R = 0. The coordinate frame 
we have arrived at is called center- oj-mass frame (CMS). 
We can also see now that 

(iri dr 

where /i — ^ — ^ is the reduced mass of the system, 

TOi -I- m2 

and hence the equation of motion can be cast in the form 



/iM. 



(2.2) 



where we have defined the radial unit vector f = -. We 

r 

can get two more conservation laws if we take the scalar 

dv 

product of Eq. ()2.2p with — and its vector product with 

dr 

r. The scalar product with — gives on the left-hand side 

dt 

dr d^r 1 d / dr 
'di 'dt^ ~ 2di \'di 



and on the right-hand side we have 



f dr 
^ ' di 



d_ 

dt 



hence 



dt\2^i r J 

where 7 = Gfj,M . This implies that the expression in 
brackets is conserved, i.e. 



— — = E = const . 

2^ r 



(2.3) 



Here the first term is the kinetic energy, the second term 
is the potential energy, and the sum of kinetic energy and 
potential energy is the total energy E, which is a constant 
of motion. Now take the cross product of Eq. (|2.2p with 
r: on the right-hand side, we get the cross product of 
coUinear vectors, which is equal to zero, hence 



d f dr\ d , . 
^dt \ dt dr 



0, 



and hence, if we define the angular momentum L by 
L = r X p , 

we get the result 

dL „ 

d^=«' 

or 

L = const , 

i.e. conservation of angular momentum. An immediate 
consequence of this conservation law is that the radius 
vector r always stays in one plane, namely the plane per- 
pendicular to L. This implies that we can without loss of 
generality choose this plane as the (xy) coordinate plane. 
The vector r is then a two-dimensional vector, 



(cos (f), sin (/)) , 



where we have defined the polar angle (p. With this nota- 
tion we can express the magnitude of angular momentum 

as 



L = r' 



dt ' 



(2.4) 



The conservation of angular momentum can be used to 
simplify the equation of motion (|2.3p . To do this we note 
that 



= (r X p)^ = r^p^ - (r • p)^ , 



hence 



T 2 

P ^^+Pr 
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where p,- = f • p is the radial component of momentum. 
Substituting into Eq. (|2.3p then gives 

2)1 r 



or, with pr — —dt, 
dr 

1 ^dr 



L2 7 



2"^ \dt J 2^r2 r 



(2.5) 



Looking back at our starting point, Eq. p.ip . we reduce 
the dimensionahty of our problem: from the simultane- 
ous differential equations of six functions of time, namely 
the six components of the position vectors ri and r2 , we 
reduce to a pair of simultaneous differential equations 
for the polar coordinates r{t) and (j){t) these equations 
contain two constants of motion, the total energy E and 
angular momentum L. Then a mass mi is moving in 
the gravitational potential $ generated by a second mass 
7712 . The vector radius and the polar angle depend on the 
time as a consequence of the star motion, i.e. r — r(t) 
and (p — (f>{t). With this choice, the velocity v of the 
mass mi can be parameterized as 

V = Vrr + VA,(h , 



The equal sign corresponds to the radial motion. For 
Emin < E < 0, the trajectory lies between a smallest 
value rmin and greatest value rmax which can be found 
from the condition E = Vef / , i.e. 



'^^\Kie) 



_L2_ 
2^' 



where the upper (lower) sign corresponds to r^ax {fmin)- 
Only for E > 0, the upper sign gives an acceptable value; 
the second root is negative and must be rejected. Let us 
now proceed in solving the differential equations (|2.4p 
and We have 



dr dr d(j) L dr L d / 1 
dt d(f> dt fir'^ d(f> ^ dcj) \r 



(2.6) 



and defining, as standard, the auxiliary variable u 
l/r, Eq. (j2.5|) takes the form 



,2 2 2fiE 
u +u - —u = — 



(2.7) 



where u = du/dcj) and we have divided by /2^. Dif- 
ferentiating with respect to 0, we get 



where the radial and tangent components of the velocity 
are, respectively. 



dr 
'dt 



We can split the kinetic energy into two terms where, due 
to the conservation of angular momentum, the second one 
is a function of r only. An effective potential energy V^f / , 



Veff = 



ill 

2^r2 
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is immediately defined. The first term corresponds to 
a repulsive force, called the angular momentum barrier. 
The second term is the gravitational attraction. The in- 
terplay between attraction and repulsion is such that the 
effective potential energy has a minimum. Indeed, differ- 
entiating with respect to r one finds that the minimum 

lies at ro = — and that 
7/i 

2 

xrrmn 

Therefore, since the radial part of kinetic energy, 

, is non-negative, the total energy must be not less than 

ymin ■ „ 



E> Er, 



'2L2 



u' (u" + u-^) =0, 



hence either u' = 0, corresponding to the circular mo- 
tion, or 



u + u = 



(2.8) 



which has the solution 



u = ^ + C cos ((f) + a) 



or, reverting the variable. 



7/i 
L2 



+ C cos [(j) + a) 



(2.9) 



which is the canonical form of conic sections in polar co- 
ordinates jsH . The constant C and a are two integration 
constants of the second order differential equation (|2.8p . 
The solution ()2.9p must satisfy the first order differential 
equation (1^ . Substituting (|^ into (P?7| we find, after 
a few algebra. 



2 _ 2/i£; /-fii 



(2.10) 



and therefore, taking account of Eq. (jll Ap . we get 
C2 > 0. This implies the four kinds of orbits given in 
Table I and in Fig. [H 



6 



c = o 


E — Emin 


circular orbits 


0< |C| < ^ 


Emin < E < 


elliptic orbits 




E = 


parabolic orbits 




E>0, 


hyperbolic orbits 



Table I. Orbits in Newtonian regime classified by the ap- 
proaching energy. 



Parabolic path 




Hypertjolic path 



and 



and since 



-7r — ' max 



7/^ r, 



max ' ' min 



— 2a, where a is the semi-major 



axis of the ellipse, one obtains 



r\<p=o 



C can be eliminated from the latter equation and Eq. 
( [^mi and then 



7 

2E 



(2.12) 



Furthermore, if we denote the distance r\^^T^/2 by I, 
the so-called semi-latus rectum or the parameter of the 
ellipse, we get 



7M 

and hence the equation of the trajectory 



(2.13) 



(2.14) 



Figure 1. Newtonian paths: in black line we have hyperbolic 
path, in blue line we have parabolic path, in red line the 
elliptical path and in ciano the circular path 



B. Circular Orbits 

Circular motion corresponds to the condition u' = 
from which one find vq — L'^ / fij where Vg// has its min- 
imum. We also note that the expression for ro together 
with Eq. (|II A|) gives 



2Er, 



(2.11) 



Thus the two bodies move in concentric circles with 
radii, inversely proportional to their masses and are al- 
ways in opposition. 



where e = 



1-; 



is the eccentricity of the ellipse. If 



we consider the major semiaxis of orbit a Eq. (|2.12p and 
the eccentric anomaly £ , the orbit can be written also as 
(see [5^ ) 

r = a{l — ecosS) , 

this equation, known as Kepler's equation, is transcen- 
dental in £, and the solution for this quantity cannot 
expressed in a finite numbers of terms. Hence, there is 
the following relation between the eccentric anomaly and 
the angle </>: 

cos£ — e 



1 — e cos £ 



D. Parabolic and Hyperbolic Orbits 



C. Elliptical Orbits 

For < |C| < hj/L'^, r remains finite for all values of 
(j). Since r{(j) + 2tt) — r{(f)), the trajectory is closed and it 
is an ellipse. If one chooses a = 0, the major axis of the 
ellipse corresponds to = 0. We get 



L2 



C 



These solutions can be dealt together. They corre- 
spond to -E > which is the condition to obtain un- 
bounded orbits. Equivalently, one has |C| > 

The trajectory is 



1(1 



(2.15) 



where e > 1. The equal sign corresponds to E' = 
. Therefore, in order to ensure positivity of r, the polar 
angle 4> has to be restricted to the range given by 
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1 + e cos (j) > . 

This means cos0 > —1, i.e. (j) ^ (~^: ^-nd the tra- 
jectory is not closed any more. For — )■ ±7r, we have 
r — > 00. The curve p.lSp . with e = 1, is a parabola. For 
e > 1, the allowed interval of polar angles is smaller than 
(j) S (— tTjTt), and the trajectory is a hyperbola. Such 
trajectories correspond to non-returning objects. Let us 
consider a semi- axis a and J- as variable, analogous to 
the elliptic eccentric anomaly £. The hyperbolic orbit is 
defined also by 

r — a(ecosh — f) , 

hence, there is the following relation between F and the 
angle <p: 

I — a(e cosh — 1) 

cos (f> = — — — . 

ea(e cosh — 1) 

Finally the parabolic orbit can be defined by the another 
relation (see [52| ) 

^ ~ 2 ' 
where P is a parameter. In this case 



COS( 



21 - P2 
P2 



As we will discuss below, this classification of orbital mo- 
tions can reveal extremely useful for the waveform sig- 
nature of gravitational radiation. Let us now take into 
account the relativistic theory of orbits. 



III. RELATIVISTIC ORBITS 

As we have seen in the above Section, the non- 
relativistic two bodies problem consists in two sub- 
problems: 

1. deriving the equations of orbital motion for two 
gravitationally interacting extended bodies, 

2. solving these equations of motion. 

In the case of widely separated objects, one can simplify 
the sub-problem by neglecting the contribution of the 
quadrupole and higher multipole momenta of the bod- 
ies to their external gravitational field, thereby approxi- 
mating the equations of orbital motion of two extended 
bodies by the equations of motion of two point masses 
located at the Newtonian center of mass of the extended 
objects. Then the sub-problem can be exactly solved as 
shown in the above Section. The two body problem in 
GR is more complicated: because of the non-linear hy- 
perbolic structure of Einstein's field equations, one is not 



sure of the good formulation of boundary conditions at 
infinity, so that the problem is not even well posed j53| . 
Moreover, since in Einstein's theory the local equations 
of motion are contained in the gravitational field equa- 
tions, it is a priori difficult to separate the problem in two 
sub-problems, as in the non-relativistic case, where one 
can compute the gravitational field as a linear functional 
of the matter distribution independently of its motion. 
Furthermore, even when one can achieve such separa- 
tion and derive some equations of orbital motion for the 
two bodies, these equations will a priori not be ordinary 
differential equations but, because of the finite velocity 
of propagation of gravity, will consist in some kind of re- 
tarded integro-differential system [s^. However, all these 
difficulties can be somehow dealt with if one resorts to 
approximation procedures and breaks the general covari- 
ance by selecting special classes of coordinates systems 

n. 

Two physically different situations, amenable to per- 
turbation treatments, have been considered in the liter- 
ature: 



1. the problem of two weakly self- gravitating, slowly 
moving, widely separated fiuid bodies which has 
been treated by the so-called FN approximation 
schemes (for references see [ssi. Is^ - lsgj ). 



2. the problem of two strongly self- gravitating, widely 
separated bodies which has been treated by match- 
ing a strong field " internal" approximation scheme 
in and near the objects to a weak field "external" 
approximations scheme outside the objects. 

The approach has been pursued both for slowly mov- 
ing objects, either BHs [60] or in general strongly self- 
gravitating objects |61|], and for strongly self-gravitating 
objects moving with arbitrary velocities [6^ |63|. In the 
latter case, equations of orbital motion were considered in 
the form of a retarded-integro-differential system which 
however could be transformed into ordinary differential 
equations and which, when attention was restricted to 
slowly moving bodies, were expanded in power series of 

— [s^, [S^l . When keeping only the first relativistic cor- 
rections to Newton's law (first post-Newtonian approxi- 
mation), it turns out that the equations of orbital mo- 
tion of widely separated, slowly moving, strongly self- 
gravitating objects depend only on two parameters (the 
Schwarzschild masses) and are identical to the equations 
of motion of weakly self-gravitating objects (when using, 
in both cases, a coordinate system which is harmonic at 
lowest order). This is, in fact, a non-trivial consequence 
of the structure of Einstein's theory 5J]. Then, in the 
next subsections, we consider the FN motion including 
secular and periodic effects at first order approximation 
and we shall show that the equations of motion can be 
written in a quasi-Newtonian form. 
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A. Relativistic Motion and conservation laws 

The relativistic case can be seen as a correction to the 
Newtonian theory of orbits [13, ■ In GR, the time is 
incorporate as a mathematical dimension, so that the 
four-dimensional rectangular perifocal coordinates are 
(x, y, z, t) and the four dimensional polar coordinates are 
(r, 0, ip, t). The first post Newtonian equations of orbital 
motion of a binary system constrain the evolution in co- 
ordinate time t of the positions ri and r2 of the two ob- 
jects. These positions represent the center of mass in the 
case of weakly self-gravitating objects (see e.g. [13, HI]) 
and the center of field in the case of strongly self gravi- 
tating objects (see (fill). They can be derived from La- 
grangian which is function of the simultaneous position 

ri(i),r2(i), and velocities Vi (i) = ^ andv2(t) — in 

at at 
a given harmonic coordinate system, and of two constant 

parameters, the Schwarzschild masses of the objects mi 

and 7712: 



the sum is over the two objects [58|, |6J]. By a Poincare 
transformation it is possible to get a PN center of mass 
frame where P pn = ^pn — 0. In this frame one has: 



JJ_ 

7711 



R 



m(™i 



7772) 



r2 



-R 



2M2c2 

/i(mi - 7772) 



7712 



dR 



2M2c2 



GM\ 



R J 



R, 



(3.4) 



where V = = vi — V2 is the istantaneous relative ve- 

dt 

locity. The problem of solving the motion of the binary 
system is then reduced to the simpler problem of solv- 
ing the relative motion in the PN center of mass frame. 
For the sake of completeness, let us write down these 
equations of motion derived from (I3.ip - p.3p . and where, 
after variation, the positions and velocities are replaced 
by their center of mass expressions p.4p : 



CpN (ri(0,r2(i),vi(i),v2(i)) = Cn + , (3.1) 



with 



and 



1 2 1 2 Gmim2 

Cn = 2"^i"i + 2"^^"2 ^ ^ — 



(3.2) 



r 1 4 

i-2 = -Tn-lUi + -7772^2 

o o 



+- 



Gmim2 
2R 



3(^1 + «2) - 7(7-1 7^2) - N^viV2 - G 



M 
(3.3) 



where we have introduced the instantaneous relative 
position vector R = ri — r2 and R = |R| while 



N 



R 

'r' 



In p.ip and (|3.2I) we used the short notations: 

vi • vi = |vi| = i;2, vi • V2 = 7'i7'2 for the ordinary Eu- 
clidean scalar products, and c is the velocity of light. The 
invariance, at the PN approximation, and modulo an ex- 
act time derivative, of CpN under spatial translations 
and Lorentz boosts implies, via Noether's theorem, the 
conservation of the total linear momentum of the system: 



dC 



PN 



dC 



PN 



PN 



9vi 9v2 ' 
and of the relativistic center of mass integral 



^-^[«MN(. + 2,-«(2Nc» 

+4{NV)Y{i' - 2) - SiNVfNv 



2NV\3iy + l))] , 
(3.5) 



where we have introduced 

/i 77ll7772 

" M 



(7771 



mass parameter 

, „ with ( < 7^ < - I . At this point 
7772)2 V ~ ~ 4,' 



it is worth to notice that in spite of the fact that it is 
in general incorrect to use, before variation, in a La- 
grangian a consequence, like Ea. (l3.4p . of the equations of 
motions, which are obtained only after variation, it turns 
out that the relative motion in PN center of mass frame, 
Eq. p.Sp . can be correctly derived from a Lagrangian 
obtained by replacing in the total Lagrangian (divided 

by fi) — >CpAr (ri, r2, vi, V2) the positions and velocities 

by their PN center of mass expressions obtained from 
(13. 4p and that moreover it is even sufficient to use the 
non-relativistic center of mass expressions: 



rijv = 


^R, 




7771 


r2jv = 


^R, 




7772 


vijv = 






7771 


V2Ar ~ 






7772 



(3.6) 



KpAT — GpN — tPpN , 

yr-^ ( 1 77ll7l2 



GpN — 



mi + -- 



1 G777l7772 

2 i?c2 



The proof goes as follows [SJ]. Let us introduce the 
following linear change spatial variables in the PN 

Lagrangian £pjv (""i ~ ''2, ^) : (ri,r2) (R,X) 



9 



With R = ri — To and X — , that is: 

M 

I"! = riAT + X, 
r2 = r2Ar + X, 

which imphes (denoting = W): 



vi = vijv 

V2 = V2JV 



w, 
w. 



Expressing 

£pN = (ri - r2, vi, V2) + £2 (ri - r2, vi, V2) , 

given by Eq. p.ip - p.3p in terms of the new variables one 
finds: 



C 



1 



PN 



1 



rC. R 



mi 



R 



W. 



V 

m2 



W 



(3.7) 



Hence one obtains as a consequence of the equations of 
the PN motion: 



_ 1 5LpN 



d_ 



d d 

dtm. 
w, 



GM 



R 



mi 



m2 



(3.8) 



1 



where in the last bracket we have discarded -MW"^ 

2. 

which gives no contribution. The first two terms in the 
[rhs) of Eq. p.Sp yield the Newtonian relative motion. 
We wish to evaluate the relativistic corrections to the rel- 
ative motion: ( — — ) ( — ^ ) in the PN center of mass 

frame. Now C2 is a polynomial in the velocities and 
therefore a polynomial in W, and from Eq. p.4p one 

sees that in the PN center of mass frame W = O ' 



Therefore as 7=- does not act on W, we see that the con- 
oR 

tributions coming from W to the rhs of Eq. (|3.8|) are of 

the second PN order O ^— ^ that we shall consistently 

neglect throughout this work. In other words one obtains 
as a consequence of the equations of the PN motion in 
the PN center of mass frame: 



_5_ 

= 



( 1 



GM 



^£2 (R,i^V,-^V 

[l& \ mi 7712 



This shows that the equations of the relative motion in 
the PN center of mass frame derive from the following 
Lagrangian: 



/:«^(R,v) = iF^ + ^ 



mi 



m2 



which happens to be obtained by replacing in the full PN 
Lagrangian, see Eq. (|3.7p above, X and W by zero, i. e. 
the original variables by Eq. (|3.6p and by dividing by /j, 
[S^ l . The explicit expression of 



Ipj^ reads: 



'-PN 



(R,V) = iy2_ 
GM 



GM 1 



R 



+ -(1-3^.)— + 



(3 + i^)V^ + i^{NV)^ - 



GM 



(3.9) 



The Lagrangian p.9|) was obtained in [gl]. The integra- 
tion of the equations (13. 5p can be done in several different 
ways. 

• A standard approach: Lagrange's method of varia- 
tion the osculating elements. 

• The Hamilton-Jacobi equation approach which, 
takes advantage of the existence of the PN La- 
grangian is the route which has been taken by Lan- 
dau and Lifshits [1^ , who worked out only the sec- 
ular precession of the periastron. 

• Another approach, based on the Maupertuis prin- 
ciple 0, which reduces the PN problem to a simple 
auxiliary Newtonian problem. 

To describe the motion, it is convenient to use the 
standard method to solve the non-relativistic two-bodies 
problem and which consists in exploiting the symmetries 
of the relative Lagrangian >Cpjv The invariance £pjv 
under time translations and space rotations implies the 
existence of four first integrals: 

V . - £5 „ and J = R X "^^^^ • 



E 



dV 



2 

GM 



dV 



GM 



3 



R X V 



(3 + v)V'^ + v{NVf - 



GM 



l + -(l-3!^)^ + (3 + i/) 



GM 
2Rc^ 



(3.10) 



(3.11) 



^ In classical mechanics, Maupertuis' principle is an integral equa- 
tion that determines the path followed by a physical system 
without specifying the time parameterization of that path. It 
is a special case of the more generally stated principle of least 
action. More precisely, it is a formulation of the equations of 
motion for a physical system not as differential equations, but as 
an integral equation, using the calculus of variations. 
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It is checked that these quantities coincide respectively 
with fi~^ times the total Noether energy and the to- 
tal Noether angular momentum of the binary system 
when computed in the PN center of mass frame [66| . 
Eq. (|3.11l) implies that the motion takes place in a 
coordinate plane, therefore one can introduce polar co- 
ordinates R — r and (/) in the plane (i.e. = 
rcos(f), Ty ~ rsin0, — 0). Then starting from 
the first integrals (|3TT0)) - (|3TT|) and using the identities: 



dt 



we obtain by iteration □ 



dr\- 
dt) 



A 



2B 

r 



C_ 



D 



d(j} _H I 
dt r'^ ' 



(3.12) 



(3.13) 



where the coefficients A, B,C, D, H, I are the following 
polynomials in E and J = |J|: 



A = 2E[l + li3.-l)^ 



E 



B = GM[ l + (7i/-6)— j , 
C = - I 1 + 2{3u - 1)4 ) + - 10) 



D = (-31^ 



MM. J 



(3.14) 



The relativistic "relative motion", i.e. the solution of 
Eq. p. 121) can be simply reduced to the integration of 
auxiliary non-relativistic radial motion. Indeed let us 
consider the following change of the radial variable: 



r = r 



D 

2C^ 



(3.15) 



where Co is the limit of C when —^0 with 

(Co = —J'^). Geometrically, the transformation which 
is expressed in polar coordinates by the equation: 
r' — r + const, (f)' — (j), is called a conchoidal trans- 
formation [54]. Taking into account the fact that D 

is C ( — ) and that we can neglect all terms of order 



O 



wc find that replacing Eq. p.l5p in Eqs. p.ip 



13.31). leads to: 



with 



C^C 



2B 

f 



BD 



(3.16) 



Then, in the case of quasi-elliptical motion {E < 0; A < 
0), r is a linear function of cosS, £ being an eccentric 

anomaly and the same is true of r = f H — - . We then 

obtain the PN radial motion in quasi-Newtonian para- 
metric form (to being a constant of integration): 



n{t — to) = £ — et sm£ . 
r = ar(l — e,- cos£) , 



(3.17) 
(3.18) 



with 



et = 



ar 



i~Af/^ 
B 



S^y 1/2 
Co )_ 



A 

B " 
1 + 



D 

'2C'o' 
AD 
2BCa 



The main difference between the relativistic radial mo- 
tion and the non-relativistic one is the appearence of two 
eccentricities: the time eccentricity et appearing in the 
Kepler equation p.l7p and the relative radial eccentric- 
ity Cr . Using (|3.14p we can express ar,er, e* and n in 
terms of E and J: 



GM 
1 + 



1, .E 



2E 
G^NP 



5 15 

2"-Y 



6) 



1/2 



et= 1 



2E 
GHP 



J2 + (-2j/ + 2) 



7 

— I 
2 

G^NP' 



17 

y 



(-2^)3/2 

GM 



1 



E 



1 + 7(^-15)3 



(3.19) 



^ In these and the following equations we neglect terms of the 
second PN order O f 



It is remarkable that a well known result of the Newto- 
nian motion is still valid at PN level: both the relative 
semi-major axis the mean motion n depend only on 
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the center of mass energy £. The same is true for the time 

27r 

of return to periastron period P ~ — . As a consequence 

n 

we can also express n in term of a^: 



GM 



(~9 + i^) 



Let us note also that the relationships between er and Ct 
are: 



1+ 31/-8K 



— = 1 



GM 



The relativistic angular motion, i.e. the solution of 
Eq. (I3.13|) can also be simply reduced to the integration 
of an auxiliary non relativistic angular motion. Let us 

first make, at O ( — ) order, the following conchoidal 



transformation: 

which transforms Eg. (13. 131) into 

dt ~ f^' 
where f can be expressed as 

f = a(l — ecos£) . 



(3.20) 



(3.21) 



Let us note also the relationship between e,- and et'. 



a = ttr 



e = 1 



/ 

' 2H ' 

^) ■ 

The differential time is given, from Eq. p. 171) by: 

dt = rt^^(l — et cos£)d£ . 

Hence we get 

iJ (1 — et cos£) 



no? (1 — ecosf )' 



-d£. 



As can be seen from Eq. p.l7|) and Eq. p.22p et and e 
differ by only small terms of order — . Now if we intro- 
duce any new eccentricity say e^ also very near et so that 

1 o , (ft + ^0) 
we can write: et = 7 r2 + e, = -— ^- e. 



(et + e^) 



with £ = O I — 1 then 



(l-etcosf)(l-e0cosf) = ( 1 - -^^^^-^ cos £ ] -e^cos^£. 



Therefore if we choose e^ such that the average of et 
and €0 is equal to e i.e. e^ = 2e — et we have 



1 



(1 — et cos£) ^ 
(1 — ecos£) 1 — e0cos£ 



O 



1 



which transforms Eq. pil Ap into a Newtonian like an- 
gular motion equation 



H 



d£ 



no? 1 — €0 cos £ 
which is easily integrated to give 

4>^4>o = KA,^{£), 



(3.23) 



being a constant of integration and where for the sake 
of simplicity we have introduced the notations: 



and 



{E) = 2 arctan 



K 



1 



H 



tan ■ 



n^il-e^i 



(3.24) 



(3.25) 



From Eq. ([5:^ and ([51^ and the definition of = 
2e — et we have: 



1 



AD 



AI 



1 



AD 



AI_\ 

BCq BH J 2BCo BH J 

then, as shown by straightforward calculations: 
G/i 



1 



2a».c^ 



2£; 
G2M2 



1 15 



)- 









J -6- 



gH'P 



and 



K 



J 



(3.26) 
(3.27) 



As it is clear from Eqs. (|3.20l) and (|3.2ip . the radial 
variable r reaches its successive minima "periastron pas- 
sages" for £ = 0, 27r,47r, ...The periastron therefore pre- 
cesses at each turn by the angle A(j) — 2'k{L — 1), which 

if J >> reduces to the well-known result j67i] : 



= 67r 



c 



dnGM 



ai?.(l 



€r 
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Contrarily to the usual approach which derives first the 
orbit by eUminating the time between Eq. (I3.12p and 
(|3.13|) before working out the motion on the orbit we 
find the orbit by eliminating S between Eq. p.lSp and 
(|3.23p - p.25p . With the aim to simplify the formulae we 
introduce the notation / for the polar angle counted from 
a periastron and corrected for the periastron precession 
i.e.: 



f 



K 



We must eliminate £ between: 

r — ar{l — Cr cosf ) , 

and 

In order to get, it is convenient to play a new conchoidal 
transformation on r writing: 



r = —ar{l — £0 cosS) + ( 1 

6(6 V £(4 



From the definition of A^^ {£) we have: 



(3.28) 



1 — 60 cos £ = 



1 



l + e^Ae^{£) 1 + e^cos/ 



Moreover we find from Eq. p.26|) that the radial dis- 
placement appearing in Eq. (j3.28p is simply 



1 ^\ - ^ 
ed> J 2c2 



so that we find the polar equation of the relative orbit 
as: 



Gfj. 
2? 



1 



"0 



1 + cos / 



2? 



This equations means that the relative orbit is the con- 
choid of a precessing ellipse, which means that it is ob- 
tained from an ellipse: r' = /(I + ecos0')~^ by a radial 
displacemnet r = r' + const together with an angular 
homothetic transformation: (f) = const ■ cj)' . Let us finally 
note that the relative orbit con also be written as: 

ar(l - el) 



1 -f Er cos /' ' 



with 



/' = / + 2(^^jsin/. 

The conservation laws and the coordinate transforma- 
tions which we have obtained here will reveal particularly 
useful to characterize the relativistic orbits, as we will see 
below. 



B. Relativistic quasi-Elliptical orbits 

The relativistic motions of each body are obtained 
by replacing the solutions for the relative motion, 
t{£),r{£),6(£), in the PN center of mass formulae Eq. 
(13. 4p (seejsj). We see first that the polar angle of the 
first object (of mass mi) is the same as the relative po- 
lar angle and that the polar angle of the second object 
(mass 7712) is simply cj) + Therefore it is sufficient 
to work out the radial motions. From Eq. (13.41) we 
have by replacing V'^ in the relativistic corrections with 

2GM 2GM GM 

+ 2E ■ 



R 



R 



m2R Gfi{mi — TO2) / R 
^ ~ M 2mc2 V ~ 

(and similar results for the second object by exchanging 
mi and 7712) which shows remarkably enough, that r can 
also be written in a quasi-Newtonian form: 



r = ar'{l — Cr' cos£) 



with 



m2 



Cr' = Er 

and where as before: 



Gmi{mi — m2 
2MaZc2 



n{t — to) = £ — et sm£ , 
^^cj)o^KA,^i£). 

The orbit in space of the first object can be written by 
using the same method as in the preceding Section for 
the relative orbit, that is: 



— Or' (1 — £0 cos £) + Gr' ( 1 — 

£0 \ . 



One finds: 



1 



Gm\ra2 
2Af2c2 



hence we find that the orbits is conchoid of a precessing 
ellipse with 



Gni\m2 
2M2c2 



1 



1 



,s(^) 



Gm\m2 
2M2c2 



Summarizing then gathering our results for the 
elliptic-like {E < 0) PN motion in the PN center of 
mass frame, we have: 



n{t — to) ~ £ — it sm£ , 



K2 arctan 



r — ar{l — Cj. cos£) , 
r' — ar>{l — e^' cos£) , 



l + £rf 



tan ■ 



(3.29) 
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with 



D. Relativistic quasi-Parabolic orbits 



GM 
~2E 



GM 



1 



1, ,E' 



1, -E' 



(3.30) 



and K, et, e^, e^, e^, a^, e^', a^' given in terms of the total 
angular momentum by unit reduced mass in the center 
of mass frame, E and J, by Eq. dXTO]) . (jXTfl) . (jX^ . 
and interchange of toi and m2 for e^'jOr'- The above 
equations are very similar to the standard Newtonian 
solutions of the non-relativistic two-body problem. 



C. Relativistic quasi-Hyperbolic orbits 

The simplest method to obtain the Post-Newtonian 
motion in the hyperbolic-like case [E > 0) consists 
simply in making, in Eas. (|3.29p - (I3.30I) . the analytic con- 
tinuation from < to i? > 0, passing below E — Q 
in the complex E plane and replacing the parameter £ 
by iJ^ (i^ = — 1). The proof that this yields to a correct 
parametric solution consists in noticing, on one hand, 
that K and the various eccentricities are analytic in i?, 
near _B = 0, and that if one replaces the parametric so- 
lution (|3.29l) - p.30l) and the corresponding expressions of 

et , 60 , Er , Br , Or , Cr' , Or' in tcrms of E and J in ( — ) and 



dt 



dt 



Eq. (13.131) are satisfied identically, modulo O — ) , and 



1 



that the resulting identities are analytic in E and £ as 
purely imaginary ones. One finds: 



n{t — to) = Ct sinh — T , 
— '/'o = K2 arctan 



1 



£0-1 
r — ar{l — Cr cosj-) , 
r' — ar'{l — Er' cosJ^) , 



, T 

tanh — 
2 



(3.31) 



where K, et, e^, Cr, Cr' are functions of E and J, as before, 
but where it has been conveninet to introduce the oppo- 
sites of analytic continuations of the semi-major axes: 



GM 
'2E 



I, ,E 



and 



miGr 

M 



and the modulus of the analytic contin- 



uation of the mean motion 
(2M) " 



GM 



1, .E 



The quasi-parabolic post-Newtonian motion {E = 0) 
can be obtained as a limit when _B — > 0. For istance, let 
us start from the quasi-elliptic solution in Eq. (|3.29|) and 
pose 



'2E 



\G^M^^ 



Taxing now the limit — >■ 0^, holding x fixed, yields 
the following parametric representation of the quasi 
parabolic motion: 



t-tn 



1 



2G^M^ 



J2 + 2 - 2i/) —X + -x^ 

3 



(3.32) 



J 



(j2_6)GiMi)i 



■2 arctan 



(J2- 6GfMi)i ' 

(3.33) 



2GM 



G^M^ 



(3.34) 



Moreover let us point out that our solutions (for the 
three cases i? < 0, £' > and £' = 0) have been 

2 [GM 

written in a suitable form when J > 6 

ever the validity of our solutions can be extended to the 



How- 



one finds that Eq. p. 121) and the square of range < 6 



by first replacing in the solutions 



— . The case of purely radial motion 
X 



for the angular motion , the second equation of (I3.29P 
and (P3T|) . and considering the Eqs. the 
function arctan by cot (at the price of a simple modifia- 
tion of 00 ) and then by making an analytic continuation 
in J. One ends up with an angular motion expressed 
by an arg coth which can also be approximatively re- 
placed by its asymptotic behaviour for large arguments: 

arg coth(X) 

( J = 0) is also obtained by taking the limit J — >■ ( at 
£,J- or respectively x fixed). Finally a parametric rep- 
resentation of the general post-Newtonian motion in an 
arbitrary (post-Newtonian harmonic) coordinate system 
is obtained from our preceding center of mass solution by 
applying a general transformation of the Poincare group 

x"" ^ Lix'' + m. 



IV. RELATIVISTIC ORBITS WITH 
GRAVITOMAGNETIC CORRECTIONS 



Using the orbital theory developed up to now for 
relativistic orbit, we have neglected terms of order 

I! 3 

— . However, we succeed in explaining, for instance, 
c 

the perihelion precession of Mercury. In cases where 
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10 < — < 10 , higher order terms hke — cannot be 

c , , c 

discarded in order to discuss consistently the problem 

of motion (see for example [68|]). In this situations, we 
are dealing with gravitomagnetic corrections. Before dis- 
cussing the theory of orbits under the gravitomagnetic ef- 
fects, let us give some insight into gravitomagnetism and 
derive the corrected metric. Theoretical and experimen- 
tal aspects of gravitomagnetism are discussed in [gl, [t^] . 

A remark is in order at this point: any theory combin- 
ing, in a consistent way, Newtonian gravity together with 
Lorentz invariance has to include a gravitomagnetic field 
generated by the mass-energy currents. This is the case, 
of course, of GR: it was shown by Lense and Thirring 
[7ll - l75| . that a rotating mass generates a gravitomag- 
netic field, which, in turn, causes a precession of plane- 
tary orbits. In the framework of the linearized weak-field 
and slow-motion approximation of GR, the gravitomag- 
netic effects are induced by the off-diagonal components 
of the space-time metric tensor which are proportional to 
the components of the matter-energy current density of 
the source. It is possible to take into account two types 
of mass-energy currents. The former is induced by the 
matter source rotation around its center of mass: it gen- 
erates the intrinsic gravitomagnetic field which is closely 
related to the angular momentum (spin) of the rotating 
body. The latter is due to the translational motion of the 
source: it is responsible of the extrinsic gravitomagnetic 
field [zl, [T^I- Let us now discuss the gravitomagnetic 
effects in order to see how they affect the orbits. 



From the Einstein field equations G^,y 
one finds 



V ^00 = , 



(4.2) 



SttG, 



hij = —^Sijp, 



16ttG , , 



(4.3) 



(4.4) 



where is the standard Laplacian operator defined 
on the fiat spacetime. To achieve Eqs. (|4.2|) - (|4.4|) . the 
harmonic condition 



has been used. 

By integrating Eqs. (I4.2I) - (|4.4I) . one obtains 



2$ 



(4.5) 



(4.6) 



A. Gravitomagnetic effects 

Starting from the Einstein field equations in the weak 
field approximation, one obtain the gravitoelectromag- 
netic equations and then the corrections on the metricQ 
(68| . The weak field approximation can be set as 



g^Luix) ^ T]^^ + hf_,^{x), \hf_,^{x)\ « 1, (4.1) 

where rj^^, is the Minkowski metric tensor and 
\hfj.uix)\ << 1 is a small deviation from it [s^. 

The stress-energy tensor for perfect - fluid matter is 
given by 



Tf"" ={p + pc^) u'^u" - pg*"" , 
which, in the weak field approximation p <^ pc^ , is 



4 



(4.7) 



The metric is determined by the gravitational Newtonian 
potential 



^x) = -G f ^^,d^x' , 
J |x-x'| 

and by the vector potential 1^', 



(4.8) 



pv 



:d^x' 



(4.9) 



given by the matter current density pv^ of the moving 
bodies. This last potential gives rise to the gravitomag- 
netic corrections. 

From Eqs. (|4.ip and (|4.5p - (14.91) . the metric tensor in 
terms of Newton and gravitomagnetic potentials is 



Notation: Latin indices run from 1 to 3, while Greek indices 
run from to 3; the flat space-time metric tensor is r?^^ = 



diag{l, —1, —1, — 1). 



2* \ 2 , 2 



2$\ 

1 Y ] ^ijdx'-dx^ . 



cdtdx-' 



(4.10) 
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From Eq. (I4.10[) it is possible to construct a variational where dl^^^^j^ = Sijdx^dx^ is the Euclidean length in- 
principle from which the geodesic equation follows. Then terval. Squaring Eq. (j4.19l) and keeping terms up to 
we can derive the orbital equations. As standard, we have order c~^, one finds 



where the dot indicates the differentiation with respect 
to the affine parameter. In order to put in evidence 
the gravitomagnetic contributions, let us explicitly cal- 
culate the Christoffel symbols at lower orders. By some 
straightforward calculations, one gets 



i Qj 



i nr 



i Oj 





1 a$ 

2 /dV^ 
\ dx^ 

1 a$ 

dx'^ 



dx^ ) 



c? \ dxi dx^ 1 



1 



9$ 



dxi 



dx^ 



dx 



(4.11) 
(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 



In the approximation which we are going to con- 
sider, we are retaining terms up to the orders $/c^ 
and /c'. It is important to point out that we are 
discarding terms like / c^)d^ /dx'^ , {V^ /e')d^/dx'' , 
{^/e')dV^/dx\ {V^/c^)dV' /dx' and of higher orders. 
Our aim is to show that, in several cases like in tight bi- 
nary stars, it is not correct to discard higher order terms 
in w/c since physically interesting effects could come out. 
The geodesic equations up to corrections are then 



dH 2 9$ dt dx^ 



da^ 
2 



dx^ da da 



'' dx^ 



dV"' \ dx* dx^ 



dx^ 



da da 



= 



(4.17) 



for the time component, and 



4$ 



C'dt' = 1 - — dlt^,i,^ + ^dx'dleucHd .(4.20) 



Inserting Eq. (|4.20p into Eq. (|4.18p . one gets, for the 
spatial components. 



d^x'' 2 a$ 



dl. 



euclid 



da^ c? dx^ \ da 
4 (dV^ . dV'^\ dl 



dx3 



dx'' 



2 9$ dx' dx'' 
c? dx' da da 

did dX'^ 



da da 







(4.21) 



Such an equation can be seen as a differential equation 
for dx'^/da which is the tangent 3- vector to the trajec- 
tory. On the other hand, Eq. (|4.2ip can be expressed in 
terms of leucUd considered as a parameter. In fact, for 
null geodesies and taking into account the lowest order 
in v/c, da is proportional to dleucUd- From Eq. (|4.17p 

/ 2 \ 
multiplied for 1 H — r-<I> , we have 



= 0. 



d f dt 2 ^dt 4 ^ ,,„,dx'\ 



da \ da 
and then 

dt 2 



da 



da J 



da 



1 



da 



(4.22) 



where, as standard, we have defined the afhne par ameter 
so that the integration constant is equal to 1 Sub- 
stituting Eq. (|4.19p into Eq. (|4.22l) . at lowest order in 
d/c, we find 



dl 



euclid 

cda 



= 1 



(4.23) 



In the weak field regime, the spatial 3-vector, tangent to 
a given trajectory, can be expressed as 



1 a$ 



dt 



da"^ (? dx^ \ da 
2 a$ dx' dx'' 4 
dx'- da da 



1 a$ , dx' dx' 

-0,;. 



dx'' da da 



dV^ ^ dV^ 

dx^ "'™ dx'' 



dt dx' p 
da da ' 



(4.18) 



for the spatial components. 

In the case of a null-geodesic, it results ds^ = da^ = 0. 
Eq. (|4.10p gives, up to the order c^^, 



4^' , , 
cdt — — —dx 



2$\ 

1 ~~ I dleuclid ■) 

J 



(4.19) 



dx'' 
da 

Through the definition 



cdx 
dip 



did 



dx'' 



dl. 



eudid 



Eq. (|4.2ip becomes 
de'' 



+ 



2 9$ 2 9$ 



l„k 



dleudid ' dx'' dx' 
4 fdV'' dV 



e e 



dx^ 



(4.24) 
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which can be expressed in a vector form as 

dp 7 A 

- [V$ - e(e • V$)] + - [e A (V A V)] 



dl 



euclid 



(4.25) 



The gravitomagnetic term is the second one in Eq. 
(|4.25|) and it is usually discarded since considered not 
relevant. This is not true if v/c is quite large as in the 
cases of tight binary systems or point masses approaching 
to black holes. 

Our task is now to achieve explicitly the trajectories, 
in particular the orbits, corrected by such effects. 



B. Gravitomagnetically corrected orbits 

Orbits with gravitomagnetic effects can be obtained 
starting from the classical Newtonian theory and then 
correcting it by successive relativistic terms. Starting 
from the above considerations (see Sec. im and IIIip we 
can see how gravitomagnetic corrections affect the prob- 
lem of orbits. Essentially, they act as a further v/c cor- 
rection leading to take into account terms up to c~^, as 
shown above. Let us start from the line element (j4.10p 
which can be written in spherical coordinates. Here 
we assume the motion of point-like bodies and then we 

can work in the simplified hypothesis $ = and 

r 

V' = It is 



ds'=^l + -^j cdt'-^l-- 

[dr^ + r^dd^ + r2 sin^ Od^''] 
8$ 

^cdt [cos 6 + sin 6 (cos 4> + sin 0)] dr 



C3 

8$ 



cdt [cos 9 (cos (t> -I- sin 0) — sin 6] rdO 



H — ^T-cdt [sin 6 (cos — sin (j))] ra 



Being, C = 1, one can multiply both members for 
1 H — — \. In the planar motion condition 9 — i: /2 , 



we obtain 



2$ 



2$ 



L2 



(cos 4> -f sin 4>)r — (cos — sin 0) 4> 



2$ 



and then, being 



2$ 



(where Rs is the 



1 



Schwarzschild radius) and it = - it is 



r 

ARsuE r ■ \ / / ■ \ 21 

(cos </? -|- sm (/3) u + (cos (fi — smip) u 

c 

^{1-Rsu) . 

By deriving such an equation, it is easy to show that, if 
the relativistic and gravitomagnetic terms are discarded, 
the Newtonian theory is recovered, being 

// , Rs 
u + u — — ^ . 

2L2 

This result probes the self-consistency of the problem. 
However, it is nothing else but a particular case since we 
have assumed the planar motion. This planarity condi- 
tion does not hold in general if gravitomagnetic correc- 
tions are taken into account. From the above Lagrangian 
(|4.26|) . it is straightforward to derive the equations of mo- 
tion 



1 



cr [rc^ + 2GM) 
-AGMt (^(cos6'(cos0-|-si 

sin 6* (cos 6 — sin 



c (rc^ + GM) I U- + sm" l r" 



cGMf^ - cGAW 



(4.27) 



As in the Newtonian and relativistic cases, from the 
line element (|4.26p . we can construct the Lagrangian 



2$\ . / 24' 

1 + — It- 1-^ 



r + r^9^+r^ sin^ 610^ 



8$ . 

- -^t [cos 9 + sin 9 (cos (j) + sin (j))] f 
8$ . 

-— t [cos 9 (cos + sin 0) — sin 9] r9 
8$ • 

- — t [sin 9 (cos (j) — sin 0)1 rcj) . 



(4.26) 



-ccote{rc^ + 2GM) Odr^ 



r2 (rc3 2GMc) 
-r (2GM CSC 6* (sin 6 - cos (j>)i + cr irc^ + GM) 



(4.28) 



1 



9 / ^ o^,^ccos6'r2 (^^^2 ^ 2GM) sir 
J.2 + 2GMc) ^ ' 

+f (AGMicos 9{cos cj) + sin cj)) - sin 9)i+ 
-2cr {rc^ + GM) o] 



(4.29) 
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corresponding to the spatial components of the 
geodesic Eq. (|4.2ip . The time component i is not nec- 
essary for the discussion of orbital motion. Being the 
Lagrangian (|4.26p £ = 1 it is easy to achieve a first in- 
tegral for f which is a natural constrain equation related 
to the energy. 

Our aim is to show how gravitomagnetic effects modify 
the orbital motion and what the parameters determining 
the stability of the problem are. As we will see, the en- 
ergy and the mass, essentially, determine the stability. 
Beside the standard periastron precession of GR, a nuta- 
tion effect is genuinely induced by gravitomagnetism and 
stability greatly depends on it. A fundamental issue for 
this study is to achieve the orbital phase space portrait. 

In Figlll the results for a given value of nutation an- 
gular velocity with a time span of 10000 steps is shown. 
It is interesting to see that, by increasing the initial nu- 
tation angular velocity, being fixed all the other initial 
conditions, we get curves with decreasing frequencies for 
f(i) and r{t). This fact is relevant to have an insight 
on the orbital motion stability (see Figl?]). The effect of 
gravitomagnetic terms are taken into account, in Fig. 
m showing the basic orbits (left) and the orbit with the 
associated velocity field in false colors (right). From a 
rapid inspection of the right panel, it is clear the sudden 
changes of velocity direction induced by the gravitomag- 
netic effects. 




Figure 2. Plots along the panellines of: r{t) (upper left), phase 
portrait of r{t) versus r{t) (bottom left), r{t) (upper right) 
and r{t) (bottom right) for a star of IMQ.The examples we 
are showing were obtained solving the system for the following 
parameters and initial conditions: /x « IMq, E = 0.95,00 ~ 
0, 00 = ^,00 = jo(^o,4'o = -Jo^o ^^'^ '^0 = - j^ and ro = 
20/i. The stiffness is evident from the trend of r{t) and r{t) 

To show the orbital velocity field, a rotation and a 
projection of the orbits along the axes of maximal energy 
can be performed. In other words, by a Singular Value 
Decomposition of the de-trended positions and velocities, 
the eigenvectors corresponding to the largest eigenvalues 
can be selected, and, of course, those representing the 
highest energy components (see Fig [3]|4l) . 

The above differential equations for the parametric or- 
bital motion are non-linear and with time-varying coef- 
ficients. In order to have a well-posed Cauchy problem, 



En= 0.95 r = 20 r = -0.1 d) = 0.01 6 = 0.001 

p p 

30 1 < < 1 < ^ 




Figure 3. Plots of basic orbits (left) The initial values are: 
A/ — IMq; En = 0.95 in mass units; ro = 20 in Schwarzschild 

radii; ^ = -^; ^ 

En= 0.95 r= 20 r = -0.1 * = 0.01 e = 0.001 

p p p 

30 I , , , , , 1 




_40l 1 1 1 1 1 1 

-30 -20 -10 10 20 30 

Figure 4. Plots of basic orbits with the associated velocity 
field. The arrows indicate the instantaneous velocities. The 
initial values are: M = IMq; E„ ~ 0.95 in mass units; 
ro = 20 in Schwarzschild radii; <^ = ^ = ^• 

we have to define: 

• the initial and final boundary condition problems; 

• the stability and the dynamical equilibrium of so- 
lutions. 

We can start by solving the Cauchy problem, as in the 
classical case, for the initial condition putting f = , 
= 0, 9 = and ^ = ^^'^ the result we get is that 
the orbit is not planar being 9^0. In this case, we are 
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En= 0.5 r =13.5 9 =0.001 

p 

0.2 I , , , , r 




-0.15' ' ' ' ' ' ' ' 

6 8 10 12 14 16 18 20 

r 

Figure 5. Breaking points examples: on the left panel, the 
first four orbits in the phase plane are shown: the red one 
is labelled 1, the green is II, the black is III and the fourth 
is IV. As it is possible to see, the orbits in the phase plane 
are not closed and they do not overlap at the orbital closure 
points; we have called this feature breaking points. In this 
dynamical situation, a small perturbation can lead the system 
to a transition to the chaos as prescribed by the Kolmogorov- 
Arnold-Moser (KAM) theorem On the right panel, it 

is shown the initial orbit with the initial (square) and final 
(circles) points marked in black. 

En= 0.5 r =13.5 9 =0.001 




r 



Figure 6. In this figure it is shown the initial orbit with the 
initial (squares) and final(circles) points marked in black. 



compelled to solve numerically the system of second or- 
der differential equations and to treat carefully the initial 
conditions, taking into account the fiigh non-linearity of 
the system. A similar discussion, but for different prob- 
lems, can be found in (79ll80l|. 

A series of numerical trials on the orbital parameters 



can be done in order to get an empirical insight on the 
orbit stability. The parameters involved in this analysis 
are the mass, the energy, the orbital radius, the initial 
values of r, 0, and the angular precession and nutation 
velocities 4> and 9 respectively. We have empirically as- 
sumed initial conditions on f, (f) and 6. 

The trials can be organized in two series, i.e. constant 
mass and energy variation and constant energy and mass 
variation. 

• In the first class of trials, we assume the mass equal 
to M — IMq and the energy i?„ (in mass units) 
varying step by step. The initial orbital radius tq 
can be changed, according to the step in energy: 
this allows to find out numerically the dynamical 
equilibrium of the orbit. We have also chosen, as 
varying parameters, the ratios of the precession an- 
gular velocity to the radial angular velocity r and 
the ratio of the nutation angular velocity 9 and the 
precession angular velocity 0. The initial condition 
on (j) has been assumed to be 0o = and the initial 
condition on 9 has been 9o ~ —■ For M = 1 (in So- 

9 1- r 

lar masses) , — = — and = , can be found out 

(j) 2 10 

two different empirical linear equations, according 
to the different values of 9, (j). One obtains a rough 
guess of the initial distance = rQ{E„) around 
which it is possible to find a guess on the equilib- 
rium of the initial radius, followed by trials and 
errors procedure. 

• In the second class of trials, we have assumed the 
variation of the initial orbital radius for different 
values of mass at a constant energy value equal to 
En = 0.95 in mass units. With this conditions, we 

f 

assume = ^^'^ assume that 9 takes the two 
values 1/2 and 1/10. One can approach the prob- 
lem also considering the mass parameterization, at 
a given fixed energy, to have an insight of the ef- 
fect of mass variation on the initial conditions. The 
masses have been varied between 0.5 and 20 Solar 
masses and the distances have been found to vary 
according to the two 3rd-order polynomial func- 
tions, according to the different values of 9 with 
respect to the mass (for details see [|68j]) 

In summary, the numerical calculations, if optimized, 
allow to put in evidence the specific contributions of grav- 
itomagnetic corrections on orbital motion. In particu- 
lar, specific contributions due to nutation and precession 
emerge when higher order terms in v/c are considered. 

The conclusion of this part of the review is that or- 
bits are highly characterized by the velocity regime of 
the moving bodies. The order of the parameter v/c de- 
termines the global shape of the trajectories. Our task is 
now to show how the motion of sources is related to the 
features of emitted GWs, in particular to their produc- 
tion and to the profile of waves. 
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Figure 7. Plots of orbits with various energy values. For 
each value of energy, four plots are shown: the first on the 
left column is the orbit, with the orbital velocity field in false 
colors. The color scale goes from blue to red in increasing 
velocity. The second on the left column is the orbit with a 
different nutation angular velocity. On the right column, the 
phase portraits r — r(r{t)) are shown. Energy varies from 0.3 
to 0.4, in mass units. The stability of the system is highly 
sensitive either to very small variation of ro or to variation on 
the initial conditions on both precession and nutation angular 
velocities: it is sufficient a variation of few percent on ro to 
induce system instability 

Part II 

Production and signature of gravitational 
waves 

The first derivation of gravitational radiation in GR is 
due to Einstein. His initial calculation [8l| was "marred 
by an error in calculation" (Einstein's words), and was 
corrected in 1918 (albeit with an overall factor of 



two error). Modulo a somewhat convoluted history (dis- 
cussed in great detail by Kennefick j83(]) owing (largely) 
to the difficulties of analyzing radiation in a nonlinear 
theory, Einstein's final result stands today as the leading- 
order " quadrupole formula" for gravitational wave emis- 
sion. This formula plays a role in gravity theory anal- 
ogous to the dipole formula for electromagnetic radia- 
tion, showing GWs arise from accelerated masses exactly 
as electromagnetic waves arise from accelerated charges. 
The quadrupole formula tells us that GWs are difficult 
to produce — very large masses moving at relativistic 
speeds are needed. This follows from the weakness of 
the gravitational interaction. A consequence of this is 
that it is extremely unlikely there will ever be an in- 
teresting laboratory source of GWs. The only objects 
massive and relativistic enough to generate detectable 
GWs are astrophysical. Indeed, experimental confirma- 
tion of the existence of GWs has come from the study 
of binary neutron star systems — the variation of the 
mass quadrupole in such systems is large enough that 
GW emission changes the system's characteristics on a 
timescale short enough to be observed. Intuitively, it is 
clear that measuring these waves must be difficult — the 
weakness of the gravitational interaction ensures that the 
response of any detector to gravitational waves is very 
small. Nonetheless, technology has brought us to the 
point where detectors are now beginning to set interest- 
ing upper limits on GWs from some sources (83 - [87| . The 
first direct detection could be, hopefully, not too far in the 
future. 



V. GRAVITATIONAL WAVES IN LINEARIZED 
GRAVITY 

The most natural starting point for any discussion of 
GWs is linearized gravity [ssl [59l . [ssj . Linearized gravity 
is an adequate approximation to GR when the spacetime 
metric, may be treated as deviating only slightly from a 
flat metric, 77^^: 

gfii^ = Vfj^i^ + ||/ip,.||<l- 

Here means "the magnitude of a typical non-zero 

component of hfj,^" . Note that the condition ^ 1 

requires both the gravitational field to be weakQ and in 
addition constrains the coordinate system to be approx- 
imately Cartesian. We will refer to /i^j^ as the metric 
perturbation; as we will see, it encapsulates GWs, but 
contains additional, non-radiative degrees of freedom as 
well. In linearized gravity, the smallness of the pertur- 
bation means that we only keep terms which are linear 
in ti^i, — higher order terms are discarded. As a con- 
sequence, indices are raised and lowered using the flat 
metric 77^,^. The metric perturbation transforms as 



^ We will work in geometrized coordinates putting c = G = 1 
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a tensor under Lorentz transformations, but not under 
general coordinate transformations. 

We now compute all the quantities which are needed to 
describe linearized gravity. The components of the afhne 
connection (Christoffel coefficients) are given by 

= -n'^" id h +d h ~ d h ) 
^ up 2 ' yyp'^pu ^ ^v'^ap ^(7"'vp} 

= \{dph^, + d,h%-d^'Kp) . 

Here means the partial derivative d/dx^. Since we 
use r]^^ to raise and lower indices, spatial indices can be 
written either in the "up" position or the "down" posi- 
tion without changing the value of a quantity: ~ fx- 
Raising or lowering a time index, by contrast, switches 
sign: /* — —ft- The Riemann tensor we construct in 
linearized theory is then given by 

(5.1) 

From this, we construct the Ricci tensor 

Rpu = R'^ppu = ^ {dpdjyhP^ + d^d^Kp - D/i^^ - d^d^h) , 

where h = /i^^ is the trace of the metric perturbation, 
and □ = dpd^ = — is the wave operator. Contract- 
ing once more, we find the curvature scalar: 

R^R\^ {dpd^h% - Uh) 

and finally build the Einstein tensor: 

= \ {dpdvh'^p + d^d^Kp - Uh^,^ - d^d^h 

-Vpudpd'^hPa + TjpuOh) . 

This expression is a bit unwieldy. Somewhat remark- 
ably, it can be cleaned up significantly by changing 
notation: rather than working with the metric per- 
turbation hfj_^, we use the trace-reversed perturbation 

hpu = hf^i, — -rj^i^h. (Notice that /z/"^ = —h, hence the 

name "trace reversed".) Replacing h^^^, with 

hpu + —rjfj^iyh in Eq. (j5.2p and expanding, we find that all 

terms with the trace h are canceled. What remains is 

Gpu = \ {daduh"^ + dPd^h^, - Dh^, - rj^^dpd^hP,) . 

(5.2) 

This expression can be simplified further by choosing an 
appropriate coordinate system, or gauge. Gauge trans- 
formations in general relativity are just coordinate trans- 
formations. A general infinitesimal coordinate transfor- 
mation can be written as x"' = x'^ + , where ^'^{x^) is 



an arbitrary infinitesimal vector field. This transforma- 
tion changes the metric via 

K"^ = hpu - 2a(^f^) , (5.3) 
so that the trace-reversed metric becomes 
h' - h' - -n h' 

A class of gauges that are commonly used in studies of 
radiation are those satisfying the Lorenz gauge condition 

d^h^, = 0. (5.4) 

(Note the close analogy to Lorentz gauge 

in electromagnetic theory, d^^A^ = 0, where is the 
potential vector.) Suppose that our metric perturbation 
is not in Lorentz gauge. What properties must satisfy 
in order to impose Lorentz gauge? Our goal is to find a 
new metric h'^^ such that — 0: 

d^'h'^, = d^hpu - d'^du^p - D^u + dudP^p 

Any metric perturbation hpi, can therefore be put into 
a Lorentz gauge by making an infinitesimal coordinate 
transformation that satisfies 

□C. = 9^ V . (5.5) 

One can always find solutions to the wave equation ()5.5p , 
thus achieving Lorentz gauge. The amount of gauge free- 
dom has now been reduced 

The amount of gauge freedom has now been reduced 
from 4 freely specifiable functions of 4 variables to 4 func- 
tions of 4 variables that satisfy the homogeneous wave 
equation D^" = 0, or, equivalently, to 8 freely specifiable 
functions of 3 variables on an initial data hypersurface. 

Applying the Lorentz gauge condition (|5.4|) to the ex- 
pression (15. 2[) for the Einstein tensor, we find that all but 
one term vanishes: 

G pu — 2^^M^ ' 

Thus, in Lorentz gauges, the Einstein tensor simply re- 
duces to the wave operator acting on the trace reversed 
metric perturbation (up to a factor —1/2). The linearized 
Einstein equation is therefore 

Ohp„ = - IGnTf^i, ; (5.6) 



^ Fairly recently, it has become widely recognized that this gauge 
was in fact invented by Ludwig Lorenz, rather than by Hendrik 
Lorentz. The inclusion of the "t" seems most likely due to confu- 
sion between the similar names; see Ref. 89| for detailed discus- 
sion. Following the practice of Griffiths ( 90], p. 42f), we bow 
to the weight of historical usage in order to avoid any possible 
confusion. 
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in vacuum, this reduces to 

nh„, 







(5.7) 



Just as in electromagnetism, the equation ()5.6|) admits a 
class of homogeneous solutions which are superpositions 
of plane waves: 

Here, oj = |k|. The complex coefficients v4^,y(k) depend 
on the wavevector k but are independent of x and t. 
They are subject to the constraint k^A^^ = (which 
follows from the Lorentz gauge condition), with k^ = 
(Wjk), but are otherwise arbitrary. These solutions are 
the gravitational waves. 






Figure 8. We show how point particles along a ring move 
as a result of the interaction with a GW propagating in the 
direction perpendicular to the plane of the ring. This figure 
refers to a wave with + polarization. 




A. Transverse traceless (TT) gauge in globally 
vacuum spacetimes 

We now specialize to globally vacuum spacetimes in 
which Tfj^y — everywhere, and which are asymptotically 
flat (for our purposes, hp_i, — > as r — > oo). Equivalently, 
we specialize to the space of homogeneous, asymptoti- 
cally flat solutions of the linearized Einstein Eq. (|5.6p . 
For such spacetimes one can, along with choosing Lorentz 
gauge, further specialize the gauge to make the metric 
perturbation be purely spatial 



'00 



Hqi = 



and traceless 



h = h/ = 0. 



(5.8) 



(5.9) 



The Lorentz gauge condition (15. 4p then implies that the 
spatial metric perturbation is transverse: 

dih^j = 0. 

This is called the transverse-traceless gauge, or TT 
gauge. A metric perturbation that has been put into 
TT gauge will be written /iJJ. Since it is traceless, there 
is no distinction between /iJJ and hJ^J . 

The conditions (|5.8p and (|5.9p comprise 5 constraints 
on the metric, while the residual gauge freedom in 
Lorentz gauge is parameterized by 4 functions that sat- 
isfy the wave equation. It is nevertheless possible to sat- 
isfy these conditions, essentially because the metric per- 
turbation satisfles the linearized vacuum Einstein equa- 
tion. When the TT gauge conditions are satisfied, the 
gauge is completely fixed. 

One might wonder why we would choose the TT gauge. 
It is certainly not necessary; however, it is extremely con- 
venient, since the TT gauge conditions completely fix all 
the local gauge freedom. The metric perturbation /iJJ 
therefore contains only physical, non-gauge information 



Figure 9. We show how point particles along a ring move 
as a result of the interaction with a GW propagating in the 
direction perpendicular to the plane of the ring. This figure 
refers to a wave with x polarization. 



about the radiation. In the TT gauge there is a close re- 
lation between the metric perturbation and the linearized 
Riemann tensor Rf^^pa- [which is invariant under the local 
gauge transformations (|5.3p by Eq. (15.11) ]. namely 



RiOjO — --ijhj^ . 

In a globally vacuum spacetime, all non-zero components 
of the Riemann tensor can be obtained from -Riojo via 
Riemann's symmetries and the Bianchi identity. In a 
more general spacetime, there will be components that 
are not related to radiation. 

Transverse traceless gauge also exhibits the fact that 
gravitational waves have two polarization components. 
For example, consider a GW which propagates in the z 
direction: hj^ = hj^{t—z) is a valid solution to the wave 
equation ^hj^ — 0. The Lorentz condition dzh^J = 
implies that hi^J{t — z) = constant. This constant 
must be zero to satisfy the condition that hab — as 
r — > cx). The only non-zero components of hj^ are then 

^JJ: h]-y:J}^i ^-nd . Symmetry and the tracefree 
condition (|5.9p further mandate that only two of these 
are independent: 



i TT 



i TT 



"-yy 



h+it-z) ; 
h^{t~z) . 



The quantities and hx are the two independent wave- 
forms of the GW (see Fig. [511) SIIH. 

To illustrate the effect of GWs on free falhng (FF) 
particles, we consider a ring of point particles initially at 
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rest with respect to a FF frame attached to the center of 
the ring, as shown in Fig. I8I9I 



h^j^ (t, z). The proper distance L between the two parti- 
cles in the presence of the GW is given by 



VI. INTERACTION OF GRAVITATIONAL 
WAVES WITH A DETECTOR 

The usual notion of "gravitational force" disappears in 
GR, replaced instead by the idea that freely falling bodies 
follow geodesies in spacetime. Given a spacetime metric 
g^u and a set of spacetime coordinates a;^, geodesic tra- 
jectories are given by the equation 



(Fx'' 



dx" dxP 
' dr dr 







where r is the proper time as measured by an observer 
travelling along the geodesic. By writing the derivatives 
in the above geodesic equation in terms of coordinate 
time t rather than proper time r, and by combining the 
fi = t (i.e coordinate) equation with the spatial, /i = j 
(i.e. spatial coordinates) equations, we obtain an equa- 
tion for the coordinate acceleration: 



d^x' 



v\Vl, + 2Tlv^ +V'>^,v^v^), (6.1) 



where v'^ = dx'^/dt is the coordinate. 

Let us now specialize to linearized theory, with the 
non-flat part of our metric dominated by a GW in TT 
gauge. Further, let us specialize to non-relativistic mo- 
tion for our test body. This implies that ^ 1, and to 
a good approximation we can neglect the velocity depen- 
dent terms in Eq. (j6.1|) : 

r^nn = . 



In linearized theory and TT gauge. 
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since Hqq 



0. Hence, we find that d^x'/dt^ = 0. 
This result could mean that the GW has no effect. This 
is not true since it just tells us that, in TT gauge, the 
coordinate location of a slowly moving, freely falling (here 
in after FF) body is unaffected by the GWs. In essence, 
the coordinates move with the waves. 

This result illustrates why, in GR, it is important to fo- 
cus upon coordinate-invariant observables (a naive inter- 
pretation of the above result would be that freely falling 
bodies are not influenced by GWs). In fact the GWs 
cause the proper separation between two FF particles to 
oscillate, even if the coordinate separation is constant. 
Consider two spatial FF particles, located at z = 0, and 
separated on the x axis by a coordinate distance Lc- Con- 
sider a GW in TT gauge that propagates down the z axis. 



dx Jl + h^J{t,z^O) 



dx 



'0 



i + -K,'it,z^o) 



(6.2) 



Notice that we use the fact that the coordinate location 
of each particle is fixed in TT gauge. In a gauge in which 
the particles move with respect to the coordinates, the 
limits of integration would have to vary. Eq. ()6.2p tells us 
that the proper separation of the two particles oscillates 
with a fractional length change 6L/L given by 



SL 1 / 



0) 



(6.3) 



Although we used TT gauge to perform this calcula- 
tion, the result is gauge independent; we will derive it 
in a different gauge momentarily. Notice that hJJ acts 
as a strain, a fractional length change. The magnitude 
/i of a wave is often referred to as the "wave strain". 
The proper distance we have calculated here is a par- 
ticularly important quantity since it directly relates to 
the accumulated phase which is measured by laser in- 
terferometric GW observatories . The "extra" phase (50 
accumulated by a photon that travels down and back the 
arm of a laser interferometer in the presence of a GW 
is 6(1) — 4tt6L I A, where A is the photon's wavelength and 
8L is the distance the end mirror moves relative to the 
beam splitteiH. We now give a different derivation of 
the fractional length change (16. 3p based on the concept 
of geodesic deviation. Consider a geodesic in spacetime 
given by x^ = z^(t\ where t is the proper time, with 
four velocity u^{t) = dz'^/dr. Suppose we have a nearby 
geodesic a;^(r) — z''{t) + ^^(t), where Lf^ir) is small. 
We can regard the coordinate displacement as a vec- 
tor L — L'^d^ on the first geodesic; this is valid to first 
order in L. Without loss of generality, we can make the 
connecting vector be purely spatial: i^'w^ = 0. Space- 
time curvature causes the separation vector to change 
with time, the geodesies will move further apart or closer 
together, with an acceleration given by the geodesic de- 
viation equation 



(6.4) 



^ This description of the phase shift only holds if L <C A, so that 
the metric perturbation does not change value very much during 
a light travel time. This condition will be violated in the high 
frequency regime for space-based GW detectors; a more careful 
analysis of the phase shift is needed in this case |92| . 
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see, e.g., Ref.[93|. This equation is valid to linear order 
in ; fractional corrections to this equation will scale as 
L/C, where C is the lengthscale over which the curvature 
varies. 

For application to GW detectors, the shortest length- 
scale £ is the wavelength A of the GWs. Thus, the 
geodesic deviation equation will have fractional correc- 
tions of order L/\. For ground-based detectors i < a 
few km, while A > SOOOkni; thus the approximation will 
be valid. For detectors with L> X (e.g. the space based 
detector LISA) the analysis here is not valid and other 
techniques must be used to analyze the detector. 

A convenient coordinate system to analyze the 
geodesic deviation equation ()6.4p is the local proper ref- 
erence frame of the observer who travels along the first 
geodesic. This coordinate system is defined by the re- 
quirements 



z\t) = 0, g,,.{t, 0) = 77,,,, ri^pit, 0) - 0, 
which imply that the metric has the form 



-dr + dx^ 



01^ 



(6.5) 



(6.6) 



Here TZ is the radius of curvature of spacetime, given by 
TZ~'^ ^ ll^pi/prrll- It also follows from the gauge con- 
ditions (|6.5I) that proper time r and coordinate time t 
coincide along the first geodesic, that u = dt and that 
Lt' = (0,L*). 

Consider now the proper distance between the two 
geodesies, which are located at .t* — and — L'^{t). 
From the metric (j6.6p we see that this proper distance 
is just |L| = \/LiLi, up to fractional corrections of order 
L^/TZ^. For a GW of amplitude h and wavelength A we 
have 7?."^ ~ h/X^, so the fractional errors are ~ hL^/X^. 
(Notice that TZ ^ C/Vh the wave's curvature scale TZ 
is much larger than the lengthscale £ characterizing its 
variations.) Since we are restricting attention to detec- 
tors with L <C A, these fractional errors are much smaller 
than the fractional distance changes ^ h caused by the 
GW. Therefore, we can simply identify |L| as the proper 
separation. 

We now evaluate the geodesic deviation equation (|6.4p 
in the local proper reference frame coordinates. From the 
conditions (16. 5p it follows that we can replace the covari- 
ant time derivative operator u'^V^ with d/{dt). Using 
u = dt and = (0, L*), we get 



d'^Ujt) 



-i?,ojo(i,0)L^(0 



(6.7) 



Note that the key quantity entering into the equation, 
RiOjOi is gauge invariant in linearized theory, so we can 
use any convenient coordinate system to evaluate it. Us- 
ing the expression (|V Ap for the Riemann tensor in terms 
of the TT gauge metric perturbation hj"^ we find 



d^L' 
~dF' 



2 dt^ 




Figure 10. Lines of force associated to the -I- (left panel) and 
X (right panel) polarizations. 



Integrating this equation using L'(t) = Lq + 5U-{t) with 
|(5L| < |Lo| gives 



5L\t) = '-hj;^{t)Ll 



(6.8) 



This equation is ideal to analyze an interferometric 
GW detector. We choose Cartesian coordinates such that 
the interferometer's two arms lie along the x and y axes, 
with the beam splitter at the origin. For concreteness, 
let us imagine that the GW propagates down the z-axis. 
Then, as discussed in Sec. lV A[ the only non-zero compo- 
nents of the metric perturbation are /iJJ = ~^yy ~ ^+ 
and h'^y = hyj = h^, where h+{t — z) and hy^{t — z) 
are the two polarization components. We take the ends 
of one of the interferometer's two arms as defining the 
two nearby geodesies; the first geodesic is defined by the 
beam splitter at x = 0, the second by the end-mirror. 
From Eq. (|6.8I) . we then find that the distances L = |L| 
of the arms end from the beam splitter vary with time as 

5L^ 



1 



L 



1 



(Here the subscripts x and y denote the two different 
arms, not the components of a vector). These distance 
variations are then measured via laser interferometry. 
Notice that the GW (which is typically a sinusoidally 
varying function) acts tidally, squeezing along one axis 
and stretching along the other. In this configuration, the 
detector is sensitive only to the -I- polarization of the GW. 
The X polarization acts similarly, except that it squeezes 
and stretches along a set of axes that are rotated with 
respect to the x and y axes by 45° . The force lines corre- 
sponding to the two different polarizations are illustrated 

in Fig. nni 

Of course, we don't expect nature to provide GWs 
that so perfectly align with our detectors. In general, 
we will need to account for the detector's antenna pat- 
tern, meaning that we will be sensitive to some weighted 
combination of the two polarizations, with the weights 
depending upon the location of a source on the sky, and 
the relative orientation of the source and the detector 

Finally, in our analysis so far of detection we have 
assumed that the only contribution to the metric per- 
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turbation is the GW contribution. However, in reality 
time-varying near zone gravitational fields produced by 
sources in the vicinity of the detector will also be present. 
From Eq. (|6.7p we see that the quantity that is actually 
measured by interferometric detectors is the space-time- 
space-time or electric-type piece Riojo of the Riemann 
tensor (or more precisely the time-varying piece of this 
within the frequency band of the detector). From the 
general expression of Riemann tensor (see [s^), we see 
that Riojo contains contributions from both hj^^ describ- 
ing GWs, and also additional terms describing the time- 
varying near zone gravitational fields. There is no way 
for the detector to separate these two contributions, and 
the time-varying near zone gravitational fields produced 
by motions of bedrock, air, human bodies, and tumble- 
weeds can all contribute to the output of the detector 
and act as sources of noise (93-l96| . 

A. The generation of gravitational waves 

GWs are generated by the matter source term on the 
right hand side of the linearized Einstein equation 

□V = -ISvrT^, , (6.9) 

cf. Eq. (|5.6p (presented here in Lorentz gauge). In this 
section we will compute the leading order contribution 
to the spatial components of the metric perturbation for 
a source whose internal motions are slow compared to 
the speed of light ( "slow- motion sources" ) . We will then 
compute the TT piece of the metric perturbation to ob- 
tain the standard quadrupole formula for the emitted 
radiation. 

Eq. (|6.9p can be solved by using a Green's function. 
A wave equation with source generically takes the form 

□/(i,x) = s(t,x) , 

where /(t,x) is the radiative field, depending on time t 
and position x, and s(i,x) is a source function. Green's 
function G'(f, x; t', x') is the field which arises due to a 
delta function source; it tells how much field is generated 
at the "field point" (t, x) per unit source at the "source 
point" (<',x'): 

□G(i, x; t', x') = S{t - t')S{x - x') . 

The field which arises from our actual source is then given 
by integrating Green's function against s(t,x): 

/(i,x) = J dt'd^x' G{t,K;t',x')s{t',x') . 

The Green's fimction associated with the wave operator 
□ is very well known (see, e.g., [S^]): 



The quantity i — |x — x'|/c is the retarded time; it takes 
into account the lag associated with the propagation of 
information from events at x to position x'. The speed 
of light c has been restored here to emphasize the causal 
nature of this Green's function; we set it back to unity 
in what follows. 

Applying this result to Eq. (16.91) . we find 

M^,x)=4/dV^Hi_^i_;^. 

Projected transverse and traceless, as already mentioned, 
the radiative degrees of freedom are contained entirely in 
the spatial part of the metric. First, consider the spatial 
part of the metric: 

We have raised indices on the right-hand side, using the 
rule that the position of spatial indices in linearized the- 
ory is irrelevant. 

We now evaluate this quantity at large distances from 
the source. This allows us to replace the factor |x — x'| 
in the denominator with r — |x|. The corresponding 
fractional errors scale as ^ L/r, where L is the size of 
the source; these errors can be neglected. We also make 
the same replacement in the time argument of Ty-: 

r,,(<-|x-x'|,x')«r.j(<-r,x'). 

Using the formula |x— x'| = r—n^x '^+0{l/r) where = 
x'^ /r, we see that the fractional errors generated by the 
replacement (jVI Ap scale as L/t, where r is the timescale 
over which the system is changing. This quantity is just 
the velocity of internal motions of the source (in units 
with c = 1), and is therefore small compared to one by 
our assumption. These replacements give 

h.,j{t, x) = ^ y" d^x' T'i{t ~ r, x') , (6.10) 

which is the first term in a multipolar expansion of the 
radiation field. 

Eq. (|6.10p almost gives us the quadrupole formula that 
describes GW emission (at leading order) . To get the re- 
maining part there, we need to manipulate this equation 
a bit. The stress-energy tensor must be conserved, which 
means c^^T'"^ = in linearized theory. Breaking this up 
into time and space components, we have the conditions 

atT™ -I- 9ir" = , 

5tr°' -I- djT' = . 
From this, it follows that 
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Multiplying both sides of this equation by x^x^ ^ we first 
manipulate the left-hand side: 

a^rOOa;^^^ = d'i {T'>°x'x^) . 



Next, manipulate the right-hand side of Eq. ()6.1ip : mul- 
tiplying by x^x^ , we obtain: 

dkdiT^^x'x^ ^ dkdi {T^^x'x^)-2dk {T'''x^ +T^^ x')+2T'^ . 

This identity is easily verified by expanding the deriva- 
tives and applying the identity diX^ — Si-' . We thus have 

92 [T^^x^x^) = dkdi {T^^x'x')-2dk {T^'^x^ +T''^x')+2T'^ 

This yields 



d X 'J- '1 7 — 



- [ d^x' {T°°x''x'') 



In going from the first to the second line, we used the 
fact that the second and third terms under the integral 
are divergences. Using Gauss's theorem, they can thus 
be recast as surface integrals; taking the surface outside 
the source, their contribution is zero. In going from the 
second to the third line, we used the fact that the in- 
tegration domain is not time dependent, so we can take 
the derivatives out of the integral. Finally, we used the 
fact that r"" is the mass density p. Defining the second 
moment Qij of the mass distribution via 



d^x' p{t,x')x"x'^ 



and combining Eqs. (|6.10p and (|6.12l) we get 



2d^Q,j{t- 



dt^ 



(6.12) 



(6.13) 



When we subtract the trace from Qij, we obtain the 
quadrupole momentum tensor: 
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To complete the derivation, we must project out the 
non-TT pieces of the right-hand side of Eq. (I6.13p . Since 
we are working to leading order in 1 /r, at each field point 
X, this operation reduces to algebraically projecting the 
tensor perpendicularly to the local direction of propaga- 
tion n = x/r, and subtracting off the trace. It is useful 
to introduce the projection tensor. 



Sij TliTlj . 



This tensor eliminates vector components parallel to n, 
leaving only transverse components. Thus, 

hfj = hkiPikPji 

is a transverse tensor. Finally, we remove the trace; what 
remains is 



lTT 



hkiPikPji — -^PijPkihki ■ 



(6.14) 



Substituting Eq. (|6.13p into (|6.14p . we obtain our final 
quadrupole formula: 



or 



2d^Qki{t-r) 
r 



dt^ 



P,k{n)P,i{n)-^-Pki{n)P,,{n) 



hJ^{t,^) = =^Qki{t-r)n,ki 



2G 
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One can now search for wave solutions of (|6.9p from 
a system of masses undergoing arbitrary motions, and 
then obtain the power radiated. The result, assuming 
the source dimensions very small with respect to the 
wavelengths, (quadrupole approximation [26]), is that 

dE . . . 

the power radiated in a solid angle f2 is 



dVl 



dE_ 
'dh 



G 



Sttc^ V dt^ 
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If one sums (I6.16P over the two allowed polarizations, 
one obtains 



E 

pol 



dE 

'dn 



1 f d^'Qu 
'2 \ dt^ 



G 

2 



di ij dj 1 

~dt^ dfi 



~ 2n. 



rik- 



1 ( d^Q. 

2 r^"^"^ 



d^Q,j d^Q 
dt^ 



dt^ 



dt^ 



-HljUk- 



d^Q 



jk 



dt^ 

(6.17) 



where h is the unit vector in the radiation direction. The 
total radiation rate is obtained by integrating (|6.17p over 
all directions of emission; the result is 



dE 
'dt 



G^QlfQ^^^'^ 
45c5 



(6.18) 



where the index (3) represents the number of differen- 
tiations with respect to time, the symbol <> indicates 
that the quantity is averaged over several wavelengths. 



B. Extension to sources with non- negligible self 
gravity 

Concerning our derivation of the quadrupole formula 
(I6.15P we assumed the validity of the linearized Einstein 
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equations. In particular, the derivation is not apphcable 
to systems with weak (Newtonian) gravity whose dynam- 
ics are dominated by self-gravity, such as binary star sys- 
tem^. This shortcoming of the above linearized-gravity 
derivation of the quadrupole formula was first pointed 
out by Eddington. However, it is very straightforward 
to extend the derivation to encompass systems with non- 
negligible self gravity. 

In full GR, we define the quantity h^'^ via 



7] 



hi"'. 



where rji^" = diag(— 1, 1, 1, 1). When gravity is weak, 
this definition coincides with our previous definition of 
/i^"^ as a trace-reversed metric perturbation. We impose 
the harmonic gauge condition 



df^iV^gn = d^hi^^ = 0. (6.19) 

In this gauge, the Einstein equation can be written as 

aa^th"" = -167r(T^'' + f"), (6.20) 

where Dflat = ?7^'^9^9y is the flat-spacetime wave oper- 
ator, and t^^" is a pseudotensor that is constructed from 
h^^'^ . Taking a coordinate divergence of this equation and 
using the gauge condition (|6.19p . stress-energy conserva- 
tion can be written 



5^(T^'^+t^'') =0. 
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Eqs. (pl^ - (p:^ and (IOT|) are precisely the same 
equations as are used in the linearized-gravity derivation 
of the quadrupole formula, except for the fact that the 
stress energy tensor T^'^ is replaced by T^"^ -I- 1^'' . There- 
fore the derivation of the last subsection carries over, with 
the modification that the formula (|6.12p for Qij is re- 
placed by 

Q^Jit) = J d?x' [r™(t,x') +<™(t,x')] 

In this equation the term t^'^ describes gravitational bind- 
ing energy, roughly speaking. For systems with weak 
gravity, this term is negligible in comparison with the 
term T^" describing the rest-masses of the bodies. There- 
fore the quadrupole formula (|6.15l) and the original defi- 
nition (|6.12p of Qij continue to apply to the more general 
situation considered here. 



C. Dimensional analysis 

The rough form of the leading GW field that we just 
derived, Eq. (|6.15p , can be deduced using simple physical 



Stress energy conservation in linearized gravity, 9''T^i/ = 0, 
forces all bodies to move on geodesies of the Minkowski metric. 



arguments. First, we define some moments of the mass 
distribution. The zeroth moment is just the mass itself: 

Mo = j pcPx = M . 

More accurately, this is the total mass-energy of the 
source. Next, we define the dipole moment: 

Ml = J pxt(fx = MLi . 

Li is a vector with the dimension of length; it describes 
the displacement of the center of mass from the origin 
we chose. As such, Mi is clearly not a very meaningful 
quantity — we can change its value simply by choosing 
a different origin. 

If our mass distribution exhibits internal motion, then 
moment of the mass current, ji = pvi, are also important. 
The first momentum is the spin angular momentum: 



Si 



pvj Xk Cijk d^x 



Finally, we look at the second momentum of the mass 
distribution: 



Mo 



p Xi Xj d^x = M Li 



where is a tensor with the dimension length squared. 

Using dimensional analysis and simple physical argu- 
ments, it is simple to see that the first moment that can 
contribute to GW emission is Af2. Consider first Mq. We 
want to combine Mq with the distance to our source, r, 
in such a way as to produce a dimensionless wavestrain 
h. The only way to do this (bearing in mind that the 
strain should fall off as l/r, and restoring factors of G 
and c) is to put 

, GMo 
h - — . 

&■ r 

Conservation of mass-energy tells us that Mq for an iso- 
lated source cannot vary dynamically. This h cannot be 
radiative; it corresponds to a Newtonian potential, rather 
than a GW. 

Let us consider now the momentum M\ . In order to get 
the right dimensions, we must take one time derivative: 

G d Ml 
h ^ . 

at r 

The extra factor of c converts the dimension of the time 
derivative to space, so that the whole expression is di- 
mensionless. Think carefully about the derivative of Mi: 



dMi 
dt 



dt 



pxid'^x= I pvid^x — Pi 



This is the total momentum of our source. Our guess for 
the form of a wave corresponding to Mi becomes 



GP 
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Also this formula cannot describe a GW. The momentum 
of an isolated source must be conserved. By boosting 
into a different Lorentz frame, we can always set P = 0. 
Terms like this can only be gauge artifacts; they do not 
correspond to radiation. Indeed, terms like (|6.22l) appear 
in the metric of a moving BH, and correspond to the 
relative velocity of the BH and the observer, jost- 

Dimensional analysis tells us that radiation from Si 
must take the form 

G d Si 
h ~ . 

at r 

Conservation of angular momentum tells us that the total 
spin of an isolated system cannot change, so we must 
reject also this term. Finally, we examine M2: 

^ G M2 
^ '^dt^~ ■ 

There is no conservation principle that allows us to reject 
this term. This is the quadrupole formula we derived 
earlier, up to numerical factors. 

In "normal" units, the prefactor of this formula turns 
out to be G/c* — a small number divided by a very 
big number. In order to generate interesting amounts of 
GWs, the variation quadrupole momentum must be enor- 
mous. The only interesting sources of GWs will be those 
which have very large masses undergoing extremely rapid 
variation; even in this case, the strain we expect from 
typical sources is tiny. The smallness of GWs reflects 
the fact that gravity is the weakest of the fundamental 
interactions. 



D. Numerical estimates 

Consider a binary star system, with stars of mass 
nil Eind 7712 in a circular orbit with separation R. The 
quadrupole moment is given by 

Q^J = (^x^xj - , (6.23) 

where x is the relative displacement, with |x| = R. We 
use the center-of-mass reference frame, and choose the 
coordinate axes so that the binary lies in the xy plane, so 
X ^ xi ~ Rcosflt, y — X2 — i? sin fit, z = X3 = 0. Let us 
further choose to evaluate the field on the z axis, so that n 
points in the z-direction. The projection operators in Eq. 
(|6.15p then simply serve to remove the zj components of 
the tensor. Bearing this in mind, the quadrupole formula 
([OS]) yields 



2Q^ 



The quadrupole moment tensor is 



cos^ Vlt — \ cos Q,t sin Vtt 



cos fit sin cos^ ili — ^ 







its second derivative is 



^2r2>i?2 



cos 2U,t 
— sin 2Vlt 




sin 2m 
-cos 2m 




The magnitude ft, of a typical non-zero component of hj!^ 



IS 



We used Kepler's 3rd la\4j to replace R with powers of 
the orbital frequency f2 and the total mass M — mi -|- 
TO2- The combination of masses here, ^iM"^/^, appears 
quite often in studies of GW emission from binaries; it 
motivates the definition of the chirp mass: 



M 



3/5 ^j2/5 _ 
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For the purpose of numerical estimate, we will take the 
members of the binary to have equal masses, so that fj, = 
M/4: 

h = . 

r 

Finally, we insert numbers corresponding to plausible 
sources: 

, ^^'^ / 1 hour \ ^/■^ / 1 kiloparsec 
2Mq) [ P J [ r 

^0.01 sec\ ^^"^ / 100 Megaparsecs 



h ~ 10"^^ 



10" 



M 
2.8 



P 



The first line corresponds roughly to the mass, distance 
and orbital period {P — 2ti/VI) expected for the many 
close binary white dwarf systems in our G alaxy. Such 
binaries are so common that they are likely to be a con- 
fusion limited source of GWs for space-based detectors, 
acting in some cases as an effective source of noise. The 
second line contains typical parameter values for binary 
neutron stars that are on the verge of spiralling together 
and merging. Such waves are targets for the ground- 
based detectors that have recently begun operations. The 
tiny magnitude of these waves illustrates why detecting 
GWs is so difficult. The emission of GWs costs energy 
and to compensate for the loss of energy, the radial sep- 
aration R between the two bodies must decrease. We 
shall now derive how the orbital frequency and GW fre- 
quency change in time, using Newtonian dynamics and 
the balance equation 



dE, 



orbit 



dt 



(6.25) 



® In units with G - 
M. 



1, and for circular orbits of radius R, R^VP ■ 
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At Newtonian order, E'orbit = —mim2/{2R). Thus, 
R = -2/3{Rn){n/n^). As long as Cl/n'^ < 1, the 
radial velocity is smaller than the tangential velocity and 
the binary's motion is well approximated by an adiabatic 
sequence of quasi-circular orbits. Eq. (j6.25p implies that 
the orbital frequency varies as 



?]2 



96 fGMn\ 



5/3 



(6.26) 



and the GW frequency /gw — 2cj, 

. _96 /3 fMy/' n/3 
/GW - y ^ [-^) /gw ■ 



(6.27) 



Introducing the time to coalescence r = icoai — t, and 
integrating Eq. (|6.27p . we get 



/gw - 130 



1.21M, 



M 



Isec 



3/8 



Hz. 



(6.28) 



where 1.21Mq is the chirp mass of a NS-NS bi- 
nary. Eq. ()6.28p predicts, e.g. coalescence times of 
~ 17min, 2sec, 1msec, for /gw ^ 10, 100, 10^ Hz. Us- 
ing the above equations, it is straightforward to compute 
the relation between the radial separation and the GW 
frequency. We find 
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M 



2.8M, 



1/3 



/lOO Hz 



gw 



2/3 
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Finally, a useful quantity is the number of GW cycles, 
defined by 



1 f*"- 1 
A/bw = - / ^(t) dt = - 

^ /«i„ ^ /fiin 



(6.30) 



Assuming rign ^ fiin, we get 



GW 
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M 



1.21M, 







-5/3 
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VII. GRAVITATIONAL WAVES FROM 
SOURCES IN NEWTONIAN MOTION 



We have now all the ingredient to characterize gravita- 
tional radiation with respect to the motion of the sources, 
i.e. with respect to different types of stellar encounters. 
Let us start with the Newtonian cases. With the above 
formalism, it is possible to estimate the amount of energy 
emitted in the form of GWs from a swtem of massive ob- 
jects interacting among them [s^ . Issj . Considering the 
quadrupole components for two bodies interacting in a 



Newtonian gravitational field, we have: 

Qxx = ^.r"^ (3 cos 20 sin - 1) , 

Qyy = /xr2(3sin20sin20 — 1) , 
Qzz = ir2/i(3cos26l-hl) , 



(7.1) 



Qxz = Qzx = r2^(|cos0sin26') , 

Qyz = Qzy = 7'^M(| siu 26* siu 0) , 

Qajy = Qij.T = '"^A* (i sin^ ^sin20) , 

where the masses mi and TO2 have polar coordinates 
{r; cos 6* cos 0, TiCosO smcj), riSmO} with i = 1,2 . We 
will work in the equatorial plane [9 = 7r/2). The origin 
of the motions is taken at the center of mass. Such com- 
ponents can be differentiated with respect to time, as in 
Eq. (|6.18p , in order to derive the amount of gravitational 
radiation in the various Newtonian orbital motions. 



A. Gravitational wave luminosity from circular and 
elliptical orbits 

The first case we are going to consider is that of closed 
circular and elliptical orbits. Using Eq. (|2.14p . let us de- 
rive the angular velocity equation 



y^G7(r7T7^i~m^(e cos (p -f 1)^ 



and then, from Eqs. ()7.ip . the third derivatives of 
quadrupolar components for the elliptical orbits are: 



XX 
d Qyy 

d^Qzz 

dt^ 

d Qxy 

dfi 



where 



Being 



= ;3(24cos0-f e(9cos20) + 11)) sin 
= -/3(24cos0-f e(13-h9cos2(/)))sin(/) 
: — 2/3esin0 

= /3(24cos(/)-^e(ll-h9cos2(/)))sin0) 
Gl{mi + 1112))^/'^ fi{e cos cj) + 1)^ 



43)q(3)«j ^ ^ [{mi+m^f ^l\l + e cos cf>)* 

(415e2 + 3(8 cos (/> + 3e cos 2(j)) 
(72 cos -I- e(70 + 27 cos 20))) sin ^^l 
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the total power radiated is given by 
dE 



dt AbcH^ 



where 



f{(j>) = {mi +TO2)^/^^(l + ecos0)'* X 

(415e^ + 3(8 cos + 3e cos 2(f>) x 
(72 cos + e(70 + 27 cos 20))) sin 0^. 

The total energy emitted in the form of gravitational ra- 
diation, during the interaction, is : 



I'OO 


dE 


Jo 


dt 



dt 



From Eq. (|2.4p . we can adopt the angle </) as a suitable 
integration variable. In this case, the energy emitted for 

01 < < 02 is 



A£;(0i,02) 



G3 



/(0) , 



and the total energy can be determined from the previous 
relation in the limits 0i — ^ and 02 ~^ i"- Thus, one has 

AE = 

where F(e) depends on the initial conditions only and it 
is given by 



Fie) 



(13824 + 1024486^ + 59412e'* + 2549e6) 
2880 



In other words, the gravitational wave luminosity strictly 
depends on the configuration and kinematics of the bi- 
nary system. 



B. Gravitational wave luminosity from parabolic 
and hyperbolic orbits 



In this case, we use Eq. (j2.15p and Eq. ()7.ip to calculate 
the quadrupolar formula for parabolic and hyperbolic or- 
bits. The angular velocity is 

= PL{ecos(j)+lf, 
and the quadrupolar derivatives are 



d^Q XX 

dt^ 

d Qyy 

~dF~ 

d^Qzz 
dt^ 

d Qxy 



where 

p^l^L^fi{ecos(f>+l)^ . 
The radiated power is given by 

dE _ _GfP_ 
dt " ^ 120c5 ^ 

[3Ue^ + (1152cos(0-h 187ecos20 

-3(80 cos 30 + 30e cos 40 + 48 cos 50 -I- 9e cos 60))e 

-192cos40-K 576], 



then 



dE Gl^L^^i^ 



7(0), 



dt 120c5 

where /(0), in this case, is 

/(0) = (314e2 + (1152cos(0-|- 187ecos20- 3(8Ocos30 
+30€ cos 40 + 48 cos 50 + 9e cos 60) )e -h 
-192cos40-H 576). 



Then using Eq. (|6.18p , the total energy emitted in the 
form of gravitational radiation, during the interaction as 
a function of 0, is given by 



Ai;(0i,02) 



1 



480c^ 



d(j) X 



= p(24 cos e(9 cos 20 11)) sin 0, 

= -p(24cos0-|-e(13-|-9cos20))sin0), 

= —2(06 sin 0, 
3 

— — -/3(ecos0 -I- l)^(5ecos0 -I- 8 cos 20 + 3ecos30), 



(Gl^L^TT (I271e^ + 24276e'' + 347686^ -t- 4608) /x^), 

and the total energy can be determined from the previ- 
ous relation in the limits 0i — >■ — tt and 02 ^ tt in the 
parabolic case. Thus, one has 

480^^^') ' 

where F(e) depends on the initial conditions only and it 
is given by 

F(e) = (l271e^ + 24276e^ + 34768e^ + 4608) . 

In the hyperbolic case, we have that the total energy is 
determined in the limits 0i — ?> ^ and 02 ^ , i.e. 

^^-"- 201600?^^')' 

where F(e) depends on the initial conditions only and is 
given by 

F(e) = [3157r(1271e^ + 24276e^ -t- 347686^ + 4608) -I- 
+16e[e[e(926704V2 - 7e(3319e2 - 32632V2e 
-t-55200)) - 383460] + 352128V2]] . 

As above, the gravitational wave luminosity strictly de- 
pends on the configuration and kinematics of the binary 
system. 
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Gravitational wave amplitude from elliptical 
orbits 



Beside luminosity, we can characterize also the GW 
amplitude starting from the motion of sources. In the 
case of a binary system and a single amplitude compo- 
nent , it is straightforward to show that 



7 11 _ 2G G(mi+m2)jj(13e cos 0+12 cos 20+£(4£-|-3 cos 30)) 



2G G(mi+™2)A'(17ecos0+12cos20+e(8e+3cos30)) 



/l22 

I 12 ^21 _ 2G G(mi +"12)^(136 cos 0+12 cos 20+e(4e+3 cos 30)) 

SO that the expected strain amplitude h ~ {h1^ + ^,22 + 
2/1^2)^^^ turns out to be 

G^(mi+TO2)^2 
" 

(3(13e cos 0+12 cos 20 + e(4e + 3 cos 30))^ 
+(17e cos 0+12 cos 20 + e(8e + 3 cos 30))^)^ , 

which, as before, strictly depends on the initial condi- 
tions of the stellar encounter. A remark is in order at 
this point. A monochromatic gravitational wave has, at 
most, two independent degrees of freedom. In fact, in the 
TT gauge, we have /i+ = hn + /122 and hx = /112 + ^21 
(see e.g. [99]). As an example, the amplitude of gravita- 
tional wave is sketched in Fig. [Tl]for a stellar encounter , 
in Newtonian motion, close to the Galactic Center. The 
adopted initial parameters are typical of a close impact 
and are assumed to be 6 = 1 AU for the impact factor 
and vq = 200 Kms~^ for the initial velocity, respectively. 
Here, we have fixed mi = TO2 = IAMq. The impact 
parameter is defined as L = bv where L is the angu- 
lar momentum and v the incoming velocity. We have 
chosen a typical velocity of a star in the galaxy and we 
are considering, essentially, compact objects with masses 
comparable to the Chandrasekhar limit (~ IAMq). This 
choice is motivated by the fact that ground-based experi- 
ments like VIRGO or LIGO expect to detect typical GW 
emissions from the dynamics of these objects or from bi- 
nary systems composed by them (see e.g. (59|). 

D. Gravitational wave amplitude from parabolic 
and hyperbolic orbits 

The single components of amplitude for a parabolic 
and hyperbolic orbits are 

^11 ^ _G^^(23ecos0+12cos20 + e(4e + 3cos30)) , 
^22 ^ Gi^L^t. (-j^yg 0+12 cos 20 + e(8e + 3 cos 30)) , 




Figure 11. The gravitational wave- forms from elliptical or- 
bits shown as function of the polar angle 0. We have fixed 
mi = m2 = IAMq. m2 is considered at rest while mi is 
moving with initial velocity vo = 200 Kms~^ and an impact 
parameter b — 1 AU. The distance of the GW source is as- 
sumed to be i? = 8 kpc and the eccentricity is e = 0.2, 0.5, 0.7. 



and then the expected strain amplitude is 



+59e^ ^ 



X (lOe^ + 9e^ cos 30 + 596^ cos 20 
(476^ + 108) ecos0 + 36)3 , 



= /i2i ^ _ 3g; L A' (4cos0 + g(cos20 + 3))sin0 



which, as before, strictly depends on the initial conditions 
of the stellar encounter. We note that the gravitational 
wave amplitude has the same analytical expression for 
both cases and differs only for the value of e which is 
e = 1 if the motion is parabolic and the polar angle range 
is G (— TT, 7r), while it is e > 1 and G (— 7r,7r) for 
hyperbolic orbits. In these cases, we have non-returning 
objects. 

The amplitude of the gravitational wave is sketched 
in Figs. [T^ and [T3] for stellar encounters close to the 
Galactic Center. As above, we consider a close impact 
and assume b = 1 AU cm, vq = 200 Kms~^ and mi — 
TO2 = 1.4M0. In summary, we can say that also in the 
case of Newtonian motion of sources, the orbital features 
characterize GW luminosities and amplitudes. 



VIII. GRAVITATIONAL WAVES FROM 
SOURCES IN RELATIVISTIC MOTION 

It is straighforward to extend the above considerations 
to orbital motions containing Post— Newtonian correc- 
tions. It is clear that GW luminosity and amplitude are 
strictly dependent on the parameter (^) considered at 
various order of approximation and, as discussed above, 
the global feature of orbits fully characterize the grav- 
itational emission. Now we study how the waveforms 
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Figure 12. The gravitational wave-forms for a parabolic 
encounter as a function of the polar angle (j). As above, 
mi = m2 = IAMq and m2 is considered at rest, mi is 
moving with initial velocity vq = 200 Kms~^ with an im- 
pact parameter 6=1 AU. The distance of the GW source is 
assumed at i? = 8 kpc. The eccentricity is e = 1. 



A. Inspiralling waveform including post-Newtonian 
corrections 

As we have shown in the above section the PN method 
involves an expansion around the Newtonian limit keet 
ing terms of higher order in the small parameter llOC 



doh 


2 






rpij 


dih 











(8.1) 



In order to be able to determine the dynamics of bi- 
nary systems with a precision acceptable for detection, 
it has been necessary to compute the force determining 
the motion of the two bodies and the amplitude of the 
gravitational radiation with a precision going beyond the 
quadrupole formula. For nonspinning BHs, the two-body 
equations of motion a nd th e GW flux are currently known 
through 3.5PN order jl03l |. Specifically if we restrict the 
discussion to circular orbits, as Eq. (|6.26p shows, there 
exists a natural adiabatic parameter fl/fl^ = 0[{v/c)^]. 
Higher -order PN corrections to Eq. (|6.26l) have been com- 
puted jl03l . Il04{ , yielding the general equation: 
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Figure 13. The gravitational wave-forms for hyperbolic en- 
counters as function of the polar angle </>. As above, we have 
fixed mi = m2 = IAAIq. m2 is considered at rest while mi is 
moving with initial velocity vq = 200 Kms~^ and an impact 
parameter 6=1 AU. The distance of the source is assumed 
at ii = 8 kpc. The eccentricity is assumed with the values 
e = 1.2,1.5,1.7 . 
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where G = 1 = c and where we define v^^ = (Mo;)^/'^ . 
The PN-order is given by ^(k/2)PN which is, up to fc = 7, 
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depend on the dynamics of binary and colliding systems 
and how relativistic corrections modulate the features of 
gravitational radiation. 
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Figure 14. We sketch the curvature potential as function of 
the tortoise coordinate r* associated to metric perturbations 
of a Schwarzschild BH. 




Figure 15. The potential peaks at the last unstable orbit for 
a massless particle (the light ring). Ingoing modes propagate 
toward the BH horizon, whereas outgoing modes propagate 
away from the source. 
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We denote L = /iX x V the Newtonian angular mo- 
mentum (with X and V, as above, the two-body center- 
of-mass radial separation and relative velocity), and 
L = L/\L\; Si — xiKiiSi and S2 = X2 m| S2 are 
the spins of the two bodies (with Si^2 unit vectors, and 
< Xi,2 < 1 for BHs) and 



S = Si + S2. 



mi 



M 



m2 



mi 



.11) 



1712 and jE = 0.577 ... is Euler's 



Finally, Sm 
constant. 

It is instructive to compute the relative contribution 
of the PN terms to the total number of GW cycles ac- 
cumulating in the frequency band of LIGO/VIRGO. In 
Table im we list the figures obtained by plugging Eq. (18. 2p 
into Eq. (|6.30|) [lOOj . As final frequency, we use the in- 
nermost stable circular orbit (ISCO) of a point particle 

44OG/(Af/M0) Hz]. 



in Schwarzschild BH [f^-^'^ 



B. The full waveform: inspiral, merger and 
ring-down 

After the two BHs merge, the system settles do wn to 
a K err BH and emits quasi- normal modes (QNMs), [l05l 
106l |. This phase is commonly known as the ring-down 
(RD) phase. Since the QNMs have complex frequencies 



totally determined by the BH's mass and spin, the RD 
waveform is a superposition of damped sinusoidals. The 
inspiral and RD waveforms can be computed analyti- 
cally. What about the merger? Since the nonlineari- 
ties dominate, the merger would be described at best and 
utterly through numerical simulations of Einstein equa- 
tions. However, before numerical relativity (NR) results 
became available, some analytic approaches were pro- 
posed. In the test-mass limit, v <^ 1, Refs. |l06l . Il07j | 
realized a long time ago that the basic physical reason 
underlying the presence of a universal merger signal was 
that when a test particle falls below 3M (the unsta- 
ble light storage ring of Schwarzschild), the GW that 
it generates is strongly filtered by the curvature poten- 
tial barrier centered around it (see Fig. [HI). For the 
equal-mass case = 1/4, Price and PuUin |108{ proposed 
the so-called close-limit approximation, which consists in 
switching from the two-body description to the one-body 
description (perturbed-BH) close to the light-ring loca- 
tion. Based on these observations, the effective-one-body 
(FOB) resummation scheme [l9| provided a first example 
of full waveform by (i) resumming the PN Hamiltonian, 
(ii) modeling the merger as a very short (instantaneous) 
phase and (iii) matching the end of the plu nge ( around 
the light-ring) with the RD phase (see Ref. |l09l | where 
similar ideas were developed also in NR). The match- 
ing was initially done using only the least damped QNM 
whose mass and spin were determined by the binary BH 
energy and angular momentum at the end of the plunge. 
An example of full waveform is given in Fig. I14I17I 

Today, with the results in NR, we are in the position 
of assessing the closeness of analytic to numerical wave- 
forms for inspiral, merger and RD phases. In Fig. [T71 
we show some first-order com paris ons between the EOB - 
analytic and NR waveforms (ll0j | (see also Ref. |llll |). 
Similar results for the inspiral phase but using PN the- 
ory [103 . 104 ] ( without r esummation) at 3.5PN order are 
given m Refs. [lis mil. So far, the agreement is qual- 
itatively good, but more accurate simulations, starting 
with the BHs farther apart, are needed to draw robust 
conclusions. 

Those comparisons are suggesting that it should be 
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Table II. Post-Newtonian contributions to the number of GW cycles accumulated from fiin — tv x 10 Hz to Slfin = ^ 
1/(6^/2 foj. binaries detectable by LIGO and VIRGO. We denote K.i ^ Si ■ L and ^ = Si ■ Sa. 
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Figure 16. GW signal from an equal-mass nonspinning BH Figure 17. GW signal from an equal-mass BH binary with 

binary as predicted at 2.5PN order by Buonanno and Damour a small spin xi = X2 = 0.06 obtained in full GR by Preto- 

(2000) in Ref. [10]. The merger is assumed almost instanta- j.jyg ^Hij^ 
neons and one QNM is included. 



possible to design purely analytic templates with the full 
numerics used to guide the patching together of the in- 
spiral and RD waveforms. This is an important avenue 
to template construction as eventually hundreds of thou- 
sands of waveform templates may be needed to extract 
the signal from the noise, an impossible demand for NR 
alone. 



IX. GRAVITATIONAL WAVES WITH 
GRAVITOMAGNETIC CORRECTIONS 

In this section, we are going to study the evolution of 
compact binary systems, formed through the capture of a 
moving (stellar) mass m by the gravitational field, whose 
source is a massive MBH of mass M where m ^ M. One 
expects that small compact objects (1 20Mq) from the 
surrounding stellar population will be captured by these 
black holes following m any-body scattering interactions 
at a relatively high rate jll2l Ill3l | . It is well known that 
the capture of stellar-mass compact objects by massive 
MBHs could constitute, potentially, a very important 



target for LISA |ll4 , Ill5j |. However, dynamics has to 
be carefully discussed in order to consider and select all 
effects coming from standard stellar mass objects inspi- 
ralling over MBHs. 

In the first part of this review, we have shown that, 
in the relativistic weak field approximation, when con- 
sidering higher order corrections to the equations of mo- 
tion, gravitomagnetic effects in the theory of orbits, can 
be particularly significant, leading also to chaotic behav- 
iors in the transient regime dividing stable from unsta- 
ble trajectories. Generally, such contributions are dis- 
carded since they are considered too small. However, in 
a more accurate analysis, this is not true and gravito- 
magnetic corrections could give peculiar characterization 
of dynamics [g^ . 

According to these effects, orbits remain rather ec- 
centric until the final plunge, and display both extreme 
rela t ivist ic perihelion precession and Lense-Thirring (tiI . 

[nuns 

precession of the orbital plane due to the spin 
of MBH, as well as orbital decay. In |ll8j |. it is illus- 
trated how the measured GW-waveforms can effectively 
map out the spacetime geometry close to the MBH. In 
[32l. [33I ■ the classical orbital motion (without relativis- 
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tic corrections in the motion of the binary system) has 
been studied in the extreme mass ratio hmit m ^ M, 
assuming the stellar system density and richness as fun- 
damental parameters. The conclusions have been that 

• the GW-waveforms have been characterized by the 
orbital motion (in particular, closed or open or- 
bits give rise to very different GW-production and 
waveform shapes); 

• in rich and dense stellar clusters, a large produc- 
tion of GWs can be expected, so that these systems 
could be very interesting for the above mentioned 
ground-based and space detectors; 

• the amplitudes of the strongest GW signals are ex- 
pected to be roughly an order of magnitude smaller 
than LISA'S instrumental noise. 

We investigate the GW emission by binary systems, in 
the extreme mass ratio limit, by the quadrupole approxi- 
mation, considering orbits affected by both nutation and 
precession effects, taking into account also gravitomag- 
netic terms in the weak field approximation of the metric. 
We will see that gravitational waves are emitted with a 
"peculiar" signature related to the orbital features: such 
a signature may be a " burst" wave-form with a maximum 
in correspondence to the periastron distance or a mod- 
ulated waveform, according to the orbit stability. Here 
we face this problem discussing in detail the dynamics 
of such a phenomenon which could greatly improve the 
statistics of possible GW sources. 

Besides, we give estimates of the distributions of these 
sources and their parameters. It is worth noticing that 
the captures occur when objects, in the dense stellar cusp 
surrounding a galactic MBH, undergo a close encounter, 
so that the trajectory becomes tight enough that orbital 
decay through emission of GWs dominates the subse- 
quent evolution. According to Refs. [zl, H^), for a typ- 
ical capture, the initial orbital eccentricity is extremely 
large (typically 1 — e ~ 10~^— 10"'^) and the initial peri- 
center distance very small (rp ^ 8 — lOOAf, where M is 
the MBH mass [ll9l|. The subsequent orbital evolution 
may (very roughly) be divided into three stages. In the 
first and longest stage the orbit is extremely eccentric, 
and GWs are emitted in short "pulses" during pericenter 
passages. These GW pulses slowly remove energy and 
angular momentum from the system, and the orbit grad- 
ually shrinks and circularizes. After ~ 10'^ — 10* years 
(depending on the two masses and the initial eccentric- 
ity) the evolution enters its second stage, where the orbit 
is sufficiently circular: the emission can be viewed as con- 
tinuous. Finally, as the object reaches the last stable or- 
bit, the adiabatic inspiral transits to a direct plunge, and 
the GW signal cuts off. Radiation reaction quickly circu- 
larizes the orbit over the inspiral phase; however, initial 
eccentricities are large enough that a substantial fraction 
of captures will maintain high eccentricity until the fi- 
nal plunge. It has been estimated [80] that about half of 
the captures will plunge with eccentricity e > 0.2. While 



individually-resolvable captures will mostly be detectable 
during the last ~ 1 — 100 yrs of the second stage (depend- 
ing on the stellar mass m and the MBH mass), radiation 
emitted during the first stage will contribute significantly 
to the confusion background. As we shall see, the above 
scenario is heavily modified since the gravitomagnetic ef- 
fects play a crucial role in modifying the orbital shapes 
that are far from being simply circular or elliptic and no 
longer closed. 

A. Gravitational waves amplitude considering 
orbits with gravitomagnetic corrections 

Direct signatures of gravitational radiation are given 
by GW-amplitudes and waveforms. In other words, the 
identification of a GW signal is strictly related to the 
accurate selection of the waveform shape by interferome- 
ters or any possible detection tool. Such an achievement 
could give information on the nature of the GW source, 
on the propagating medium, and, in principle, on the 
gravitational theory producing such a radiation jl20l | . 

Considering the formulas of previous Section, the GW- 
amplitude can be evaluated by 

R being the distance between the source and the observer 
and, due to the above polarizations, {j, fc} — 1, 2. 

From Eq. (|9.1l) . it is straightforward to show that, for 
a binary system where m <^ M and orbits have gravito- 
magnetic corrections, the Cartesian components of GW- 
amplitude are 

/i^^ = 2/i[(3 cos^ sin^ 6 - l)r^ + 6r{e cos^ sin 29 
— (/)sin^ 9 sin 2(f>)r + r{{3 cos^ i^sin^ 9 — l)f 
+3r{9'^ cos 29 cos^ (j) ~ 4^9 sin 29 sin 2(1) 
— sin 9 {sin 9 {(j)^ cos2(/) -I- 0cos(/)sin^) — 9 cos6'cos^ '/')))] i 

= 2fi[{3sin^9sin^(j)- l)r'^ + 6r{^sin2(j)sin^ 9 
+9 sin 29 sin^ (j))r + +r((3 sin^ 9 sin^ (jj — l)r 
+3r{9^ cos 29 sin^ (j) + (j)9 sin 29 sin 20 
-I- sin6'(^ cos6'sin^ (/> -I- sin9{4'^ cos2(/) + cj) cos (fi sin </>))))], 

^xy ^ j^yx ^ 3^[cos 20 sin 61(4^0 COS 61 -h 4>sin9y 
+2r (20 cos 20 sin^ 6* -I- ^ sin 29 sin 20)r 
+ ^ sin 20(2f sin^ 9 + r{29^ cos 29 - 40^ sin^ 9 
+9 sin 29))r + sin'- 9 sin 20] , 

where we are assuming geometrized units. The above 
formulas have been obtained from Eas. (l4.27l) . (|4.28l) . 
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(|4.29[) The gravitomagnetic corrections give rise to signa- 
tures on the GW-ampHtudes that, in the standard New- 
tonian orbital motion, are not present (see for example 
[3^ [33I). On the other hand, as discussed in the Intro- 
duction, such corrections cannot be discarded in peculiar 
situations as dense stellar clusters or in the vicinity of 
galaxy central regions. We are going to evaluate these 
quantities and results are shown in Figs. [THJ [111 HOI HIl 



B. Numerical results 

Now we have all the ingredients to estimate the effects 
of gravitomagnetic corrections on the GW-radiation. 
Calculations have been performed in geometrized units 
in order to evaluate better the relative corrections in ab- 
sence of gravitomagnetism. For the numerical simula- 
tions, we have assumed the fiducial systems constituted 
by a TO = 1.4Mq neutron star or to — IOA/q massive 
stellar object orbiting around a MBH Af ~ 3 x lO^M© 
as SgrA*. In the extreme mass-ratio limit, this means 

mM 

that we can consider /i = of about n w \AMq 

m + M 

and ~ IOMq. The computations have been performed 
starting with orbital radii measured in the mass unit 
and scaling the distance according to the values shown 
in Table I. As it is possible to see in Table I, start- 
ing from ro = 20^ up to 2500/i, the orbital eccentric- 

ity e = evolves towards a circular orbit. In 

Table I, the GW-frequencies, in mHz, as well as the h 
amplitude strains and the two polarizations /i+ and hx 
are shown. The values are the mean values of the GW 
h amplitude strains (h = ) and the maxima of 

the polarization waves (see Figs. [20l and [21]) . In Fig. [22| 
the fiducial LISA sensitivity curve is shown 13] consid- 
ering the confusion noise produced by White Dwarf bi- 
naries (blue curve). We show also the h amplitudes (red 
diamond and green circles for fi w IAMq and w IOMq 
respectively). It is worth noticing that, due to very high 
Signal to Noise Ratio, the binary systems which we are 
considering result extremely interesting, in terms of prob- 
ability detection, for the LISA interferometer (see Fig. 



X. RATE AND EVENT NUMBER 
ESTIMATIONS IN DENSE STELLAR SYSTEMS 

At this point, it is important to give some estimates 
of the number of events where gravitomagnetic effects 
could be a signature for orbital motion and gravitational 
radiation. From the GW emission point of view, close 
orbital encounters, collisions and tidal interactions have 
to be dealt on average if we want to investigate the gravi- 
tational radiation in a dense stellar system. On the other 
hand, dense ste llar r egions are the favored target for LISA 
interferometer (ll5{ so it is extremely useful to provide 



suitable numbers before its launching. 

To this end, it is worth giving the stellar encounter rate 
producing GWs in astrophysical systems like dense glob- 
ular clusters or the Galactic Center. In general, stars are 
approximated as point masses. However, in dense regions 
of stellar systems, a star can pass so close to another that 
they raise tidal forces which dissipate their relative or- 
bital kinetic energy and the Newtonian mechanics or the 
weak field limit of GR cannot be adopted as good ap- 
proximations. In some cases, the loss of energy can be so 
large that stars form binary (the situation which we have 
considered here) or multiple systems; in other cases, the 
stars collide and coalesce into a single star; finally stars 
can exchange gravitational interaction in non-returning 
encounters. 

To investigate and parameterize all these effects, one 
has to compute the collision time tcou, where 1/tcoU is 
the collision rate, that is, the average number of physical 
collisions that a given star suffers per unit time. As a 
rough approximation, one can restrict to stellar clusters 
in which all stars have the same mass m. 

Let us consider an encounter with initial relative veloc- 
ity vq and impact parameter b. The angular momentum 
per unit mass of the reduced particle is L = bvQ. At the 
distance of closest approach, which we denote by rcou, 
the radial velocity must be zero, and hence the angular 
momentum is L — rcoiiVmax, where Vmax is the relative 
speed at rcoii- It is easy to show that [1^ 



AGmrr 



coll 



(10.1) 



If we set VcoH equal to the sum of the radii of the two 
stars, a collision will occur if the impact parameter is 
less than the value of b, as determined by Eq. (|10.ip . 

The function f{'Va)(P'Va gives the number of stars per 
unit volume with velocities in the range + dP^a- The 
number of encounters per unit time with impact param- 
eter less than &, which are suffered by a given star, is 
/(va)(i^Va times the volume of the annulus with radius 
b and length vq, that is. 



f{Va)nb'^Vod^Va, 



(10.2) 



where vq — |v — Va| and v is the velocity of the consid- 
ered star. The quantity in Eg. (110.21) is equal to l/tcoii for 
a star with velocity v: to obtain the mean value of l/tcou, 
we average over v by multiplying (|10.2|) by /(v) /i>, where 
V — ^ f{v)(P\r is the number density of stars and the in- 
tegration is over d^v. Thus 

1 

tcoU 



AGmrcoii 



TcoU |v - + ) d^vdSv, .(10.3) 



Replacing the variable by V = v — V(j, the argument 



of the exponential is then — 



1^2 



/a^, and 
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Orbital parameters ot the basic orbit: T=11317 sec, e=0.64, fGW=0-2 mHz 
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Figure 18. Plots of zi^o{t) {left upper panel) and ZQrav{t) (right upper panel). It is interesting to see the differences of about five orders 
of magnitude between the two plots. At the beginning, the effect is very small but, orbit by orbit, it grows and, for a suitable interval of 
coordinated time, the effect cannot be neglected (see the left bottom panel in which the differences in x and y, starting from the initial 
orbits up to the last ones, by steps of about 1500 orbits, are reported). The internal red circle represents the beginning, the middle one is 
the intermediate situation (green) and the blue one is the final result of the correlation between Ax versus Ay, being Ax = xcrav — ^NO 
and Ay = ycrav ~ UNO- the bottom right, it is shown the basic orbit. 



Table III. GW-amplitudes and frequencies as function of eccentricity e, reduced mass /i, orbital radius ro for the two cases of 
fiducial stellar objects m ~ IAMq and m ~ IOMq orbiting around a MBH of mass M ~ 3 X 10® M©. 
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if we replace the variable v by Vcm = v — -V (the center 
of mass velocity), then one has 

tcoii stt^o-*' J \ y J 

(10.4) 

The integral over Vcm is given by 

j e-^'^/'^^d^v^ = ^'^a^ . (10.5) 

Thus 
tcoii 2o-^ 

(10.6) 



The integrals can be easily calculated and then we find 

^ A^/nvar^^ii + . (10.7) 

tcoll cr 

The first term of this result can be derived from the ki- 
netic theory. The rate of interaction is z^S {V), where E 
is the cross-section and {V) is the mean relative speed. 
Substituting E — 7rr^^;; and (V) = Aa/y/n (which is ap- 
propriate for a Maxwellian distribution with dispersion 
a) we recover the first term of (|10.7p . The second term 
represents the enhancement in the collision rate by grav- 
itational focusing, that is, the deflection of trajectories 
by the gravitational attraction of the two stars. 

If r* is the stellar radius, we may set VcoU = 2r*. It 
is convenient to introduce the escape speed from stellar 
/ 2Gm 

surface, = W , and to rewrite Ea. (|10.7p as 
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Figure 19. Plot of the differences of total gravitational wave- 
form h, with and without the gravitomagnetic orbital correc- 
tion for a neutron star of 1.4Mq orbiting around a MBH . 
The waveform has been computed at the Earth-distance from 
SgrA* (the central Galactic Black Hole). The example we 
are showing has been obtained solving the systems for the 
following parameters and initial conditions: /x ~ 1.4Mq, ro, 
E = 0.95, 00 = 0, 00 = f , ^0 = 0, </)o = -j^ro and fo = "t^- 
It is worth noticing that frequency modulation gives cumula- 
tive effects after suitable long times. 
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tcoU 



4^2 ) 



where 



e 



Gm 
2^ 



(10.8) 



(10.9) 



is the Safronov number j25[. In evaluating the rate, we 
are considering only those encounters producing gravita- 
tional waves, for example, in the LISA range, i.e. between 
lO"** and 10~^ Hz (see e.g. [12l| ). Numerically, we have 



r ~ 5.5 X 10 



e << 1 



-10 



lOkms- 



KUA^J \ R 



/lGpc\ , 
' yrs ^ 



(10.10) 



— -^"(lO^)^(lOh^)(^) 

(^)'yrs- e»i (1 



(10.11) 



If 9 >> 1, the energy dissipated exceeds the relative 
kinetic energy of the colliding stars, and the stars coalesce 
into a single star. This new star may, in turn, collide and 
merge with other stars, thereby becoming very massive. 
As its mass increases, the collision time is shorten and 
then there may be runaway coalescence leading to the 
formation of a few supermassive objects per clusters. If 
<< 1, much of the mass in the colliding stars may be 
liberated and forming new stars or a single supermassive 
objects (see jl22l . Il23| ). Both cases are interesting for 
LISA purposes. 



Note that when one has the effects of quasi-collisions 
(where gravitomagnetic effects, in principle, cannot be 
discarded) in an encounter of two stars in which the min- 
imal separation is several stellar radii, violent tides will 
raise on the surface of each star. The energy that excites 
the tides comes from the relative kinetic energy of the 
stars. This effect is important for 9 >> 1 since the loss 
of small amount of kinetic energy may leave the two stars 
with negative total energy, that is, as a bounded binary 
system. Successive peri-center passages will dissipates 
more energy by GW radiation, until the binary orbit is 
nearly circular with a negligible or null GW radiation 
emission. 

Let us apply these considerations to the Galactic Cen- 
ter which can be modelled as a system of several compact 
stellar clusters, some of them similar to very com pact 
globular clusters with high emission in X-rays jl24| . 

For a typical globular cluster around the Galactic Cen- 
ter, the expected event rate is of the order of 2 x 10~^ 
yrs~^ which may be increased at least by a factor ~ 100 
if one considers the number of globular clusters in the 
whole Galaxy. If the stellar cluster at the Galactic Cen- 
ter is taken into account and assuming the total mass 
M ~ 3 X 10^ M0, the velocity dispersion a ~ 150 km s~^ 
and the radius of the object i? ~ 10 pc (where 9 = 4.3), 
one expects to have ~ 10~^ open orbit encounters per 
year. On the other hand, if a cluster with total mass 
M ~ 10*5 Mq, cr ~ 150 km s"^ and i? ~ 0.1 pc is consid- 
ered, an event rate number of the order of unity per year 
is obtained. These values could be realistically achieved 
by data coming from the forthcoming space interferom- 
eter LISA. As a secondary effect, the above wave-forms 
could constitute the "signature" to classify the different 
stellar encounters thanks to the differences of the shapes 
(see Figs. [^and l^ . 



XI. DISCUSSION, CONCLUSIONS AND 
PERSPECTIVES 



We have considered the two-body problem in Newto- 
nian and relativistic theory of orbits in view of char- 
acterizing the gravitational radiation, starting from the 
motion of the sources. We have reported several re- 
sults concerning the equations of motion, and the asso- 
ciated Lagrangian formulation, of compact binary sys- 
tems. These equations are necessary when construct- 
ing the theoretical templates for searching and analyz- 
ing the GW signals from inspiralling compact binaries in 
VIRGO-LIGO and LISA type experiments. Considering 
the two-body problem, we mean the problem of the dy- 
namics of two structureless, non-spinning point-particles, 
characterized by solely two mass parameters mi and m2, 
moving under their mutual, purely gravitational interac- 
tion. Surely this problem, because of its conceptual sim- 
plicity, is among the most interesting ones to be solved 
within any theory of gravity. Actually, there are two as- 
pects of the problem: the first sub-problem consists into 
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Figure 20. Plots along the panel lines from left to right of field velocities along the axes of maximum covariances, total 
gravitational emission waveform h and gravitational waveform polarizations h+ and hx for a neutron star of 1.4Mq. The 
waveform has been computed for the Earth-distance from Sagittarius A (the central Galactic Black Hole SgrA*). The plots 
we are showing have been obtained solving the system for the following parameters and initial conditions: ^ ~ IAMq, 
E = 0.95, 4>o = 0, 9o = |-, do = 0, </>o = —Jq'T'q and fo ~ ^y^- From top to bottom of the panels, the orbital radius is 
ro = 20m, 1500/x, 2500^. See also Table I. 



obtaining the equation of the binary motion, the second 
is to find the (hopefully exact) solution of that equation. 
We referred to the equation of motion as the explicit ex- 
pression of the acceleration of each of the particles in 
terms of their positions and velocities. It is well known 
that in Newtonian gravity, the first of these sub-problems 
is trivial, as one can easily write down the equation of 
motion for a system of N particles, while the second one 
is difficult, except in the two-body case N = 2, which 
represents, in fact, the only situation amenable to an ex- 
act treatment of the solution. In GR, even writing down 
the equations of motion in the simplest case = 2 is 
difficult. Unlike in Newton's theory, it is impossible to 
express the acceleration by means of the positions and ve- 
locities, in a way which would be valid within the exact 
theory. Therefore we are obliged to resort to approx- 
imation methods. Let us feel reassured that plaguing 
the exact theory of GR with approximation methods is 



not a shame. It is fair to say that many of the great 
successes of this theory, when confronted to experiments 
and observations, have been obtained thanks to approx- 
imation methods. Furthermore, the beautiful internal 
wheels of GR also show up when using approximation 
methods, which often deserve some theoretical interest 
in their own, as they require interesting mathematical 
techniques. Here we have investigated the equation of the 
binary motion in the post-Newtonian approximation, i. e. 
as a formal expansion when the velocity of light c tends 
to infinity. As a consequence of the equivalence princi- 
ple, which is incorporated by hand in Newton's theory 
and constitutes the fundamental basis of GR, the accel- 
eration oi particlel should not depend on mi (nor on its 
internal structure) , in the test-mass limit where the mass 
mi is much smaller than m2. This is, of course, satisfied 
by the Newtonian acceleration, which is independent of 
mi, but this leaves the possibility that the acceleration 
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Figure 21. Plots along the panel lines from left to right of field velocities along the axes of maximum covariances, total 
gravitational emission waveform h and gravitational waveform polarizations h+ and hx for a Black Hole (BH) of IOMq. The 
waveform has been computed for the Earth-distance to SgrA*. The plots we are showing have been obtained solving the system 
for the following parameters and initial conditions: fi ~ IOMq, E = O.95,0o = 0, So = f ,^0 ~ 0,(/!)o = "iij'^o and ro = 
From top to bottom of the panels, the orbital radius is ro = 20/i, 1000/^, 2500/^. See also Table I 



of the particlel, in higher approximations, does depend 
on mi, via the so-called self- forces, which vanish in the 
test-mass limit. Indeed, this is what happens in the post- 
Newtonian and gravitomagnetic corrections, which show 
explicitly many terms proportional to (powers of) mi. 
Though the approximations and corrections to the or- 
bits are really a consequence of GR, they should be in- 
terpreted using the common-sense language of Newton. 
That is, having chosen a convenient general-relativistic 
(Cartesian) coordinate system, like the harmonic coordi- 
nate system adopted above, we have express the results 
in terms of the coordinate positions, velocities and ac- 
celerations of the bodies. Then, the trajectories of the 
particles can be viewed as taking place in the absolute 
Euclidean space of Newton, and their (coordinate) ve- 
locities as being defined with respect to absolute time. 
Not only this interpretation is the most satisfactory one 
from a conceptual point of view, but it represents also 
the most convenient path for comparing the theoretical 
predictions and the observations. For instance, the So- 



lar System dynamics at the first post-Newtonian level is 
defined, following a recent resolution of the International 
Astronomical Union, in a harmonic coordinate system, 
the Geocentric Reference System (GRS), with respect to 
which one considers the absolute motion of the planets 
and satellites. But because the equations come from GR, 
they are endowed with the following properties, which 
make them truly relativistic. 

• The one-body problem in GR corresponds to the 
Schwarzschild solution, so the equations possess 
the correct perturbative limit, that given by the 
geodesies of the Schwarzschild metric, when the 
mass of one of the bodies tends to zero. 

• Because GR admits the Poincare group as a global 
symmetry (in the case of asymptotically flat space- 
times), the harmonic-coordinate equations of mo- 
tion stay invariant when we perform a global 
Lorentz transformation 

• Since the particles emit gravitational radiation 
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Figure 22. Plot of estimated mean values of GW-emission in 
terms of strain h for two binary sources at the Galactic Center 
SgrA* with reduced mass fj, ~ IAMq (red diamonds) and 
fi ~ IOA/0 (green circles). The blue line is the foreseen LISA 
sensitivity curve. The waveforms have been computed for 
the Earth-distance to SgrA*. The examples we are showing 
have been obtained solving the systems for the parameters 
and initial conditions reported in Figs. 1201 1211 and in Table I. 



there are some terms in the equations which are 
associated with radiation reaction. These terms ap- 
pear at the order 2.5PN or c~5 that we discarded 
in our discussion (where 5 = 2s-|-l,s = 2 being 
the heUcity of the graviton). They correspond to 
an odd- order PN correction, which does not stay 
invariant in a time reversal. By contrast, the even- 
orders, as IPN, correspond to a dynamics which is 
conservative. 

• GR is a non-linear theory (even in vacuum), and 
some part of the gravitational radiation which was 
emitted by the particles in the past scatters off 
the static gravitational field generated by the rest- 
masses of the particles, or interacts gravitationally 
with itself. 

From all these considerations, the post-Newtonian equa- 
tions were also obtained, for the motion of the centers 
of mass of extended bodies, using a technique that can 
be qualified as more physical than the surface-integral 
method, as it takes explicitly into account the structure 
of the bodies. 

Particularly interesting is considering gravitomagnetic 
effects in the geodesic motion. In particular, one can 
consider the orbital effects of higher-order terms in v/c 
which is the main difference with respect to the stan- 
dard approach to the gravitomagnetism. Such terms are 
often discarded but, as we have shown, they could give 
rise to interesting phenomena in tight binding systems 
as binary systems of evolved objects (neutron stars or 



black holes). They could be important for objects falling 
toward extremely massive black holes as those seated in 
the galactic centers [zl, [sO] ■ The leading parameter for 
such correction is the ratio v/c which, in several physical 
cases cannot be simply discarded . For a detailed dis- 
cussion see for example |l25l - ll28| . A part the standard 
periastron precession effects, such terms induce nutations 
and are capable of affecting the stability basin of the or- 
bital phase space. As shown, the global structure of such 
a basin is extremely sensitive to the initial angular ve- 
locities, the initial energy and mass conditions which can 
determine possible transitions to chaotic behaviors. De- 
tailed studies on the transition to chaos could greatly aid 
in gravitational wave detections in order to determine the 
shape, the spectrum and the intensity of the waves (for 
a discussion see (l29l Il30l |). 

In the second part of this review, we have summarized 
many of the most important topics in the theory of GWs. 
Linearized theory as described in is adequate to describe 
the propagation of GWs and to model the interaction of 
GWs with our detectors. A variety of formalisms have 
been developed. 

• Newtonian theory The emission of gravitational 
waves from stellar encounters in Newtonian regime 
interacting on elliptical, hyperbolic and parabolic 
orbits is studied in the quadrupole approxima- 
tion. Analytical expressions are then derived for 
the gravitational wave luminosity, the total energy 
output and gravitational radiation amplitude pro- 
duced in tight impacts where two massive objects 
closely interact at an impact distance of lAU. 

• Post- Newtonian theory. PN theory is one of the 
most important of these formalisms, particularly 
for modeling binary systems. Roughly speaking, 
PN theory analyzes sources using an iterated ex- 
pansion in two variables: The "gravitational po- 
tential", $ ^ M/r, where M is a mass scale and r 
characterizes the distance from the source; and ve- 
locities of internal motion, v. (In linearized theory, 
we assume $ is small but place no constraints on v.) 
Newtonian gravity emerges as the first term in the 
expansion, and higher order corrections are found 
as the expansion is iterated to ever higher order. 
Our derivation of the quadrupole formula gives the 
leading order term i n the PN expansion of the emit- 
ted radiation. See |131| and references therein for 
a comprehensive introduction to and explication of 
this subject. 

• Gravitomagnetic corrections. The gravitomagnetic 
effect could give rise to interesting phenomena in 
tight binding systems such as binaries of evolved 
objects (NSs or BHs). The effects reveal particu- 
larly interesting if — is in the range (10~^ 10~^)c. 

c 

They could be important for objects captured and 
falling toward extremely massive black holes such 
as those at the Galactic Center. Gravitomagnetic 
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orbital corrections, after long integration time, in- 
duce precession and nutation and then modifica- 
tion on the wave-form. In principle, GW emission 
could present signatures of gravitomagnetic correc- 
tions after suitable integration times in particular 
for the on going LISA space laser interferometric 
GW antenna. 

To conclude, Henri Poincare (l32| once remarked that 
real problems can never be classified as solved or unsolved 
ones, but that they are always more and less solved. This 
remark applies particularly well to the problem of mo- 
tion, which has had chequered history. Even the New- 
tonian problem of motion, which appeared to well un- 
derstood after the develo pmen t of the powerful meth- 
ods of classical mechanics [l33j | embarked on an entirely 
new career after work of Po inca r e wh ich has led to many 
further developments (see |134| . Il35{ ). The Einsteinian 
problem of motion has not even reached a classical stage 
where the basic problems appear as well understood. At 
first sight the best developed approximation method in 
GR, the PN one, would seem to constitute such classical 
stage, but the literature on the PN problem of motion 
is full of repetitions, errors or ambiguities. We was to 



conclude this review by giving a list of issues that need 
to be clarified. We renounced this project because, if one 
wishes to look at the work done with a critical eye, nearly 
all aspects of the problem of motion and GWs need to 
be thoroughly re-investigates for mathematical, physical 
or conceptual reasons; so that the list of open problems 
would, consistent with the remark of Poincare. One thing 
is certain: the problem of motion and GWs is no longer a 
purely theoretical problem, tanks to the improvement in 
the precision of positions measurements in the solar sys- 
tems, and to the discovery of the binary pulsar 1913-1-16 
which is a relativistic laboratory; the problem has become 
an important tool of modern astrophysics. It is therefore 
of some urgency, not only to complete and unify the work 
already done, but also to develop new approaches in or- 
der to aim both formal and conceptual clarification of the 
basic issues, and to obtain more accurate explicit results. 
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